arXiv:l506.01402v2 [astro-ph.EP] 10 Sep 2015 


Version September 11, 2015 

Preprint typeset using I^-TjeX style emulateapj v. 08/22/09 


WAS COMET C/1945 XI (DU TOIT) A DWARF, SOHO-LIKE KREUTZ SUNGRAZER? 

Zdenek Sekanina 1 & Rainer Kracht 2 

1 -Jel Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, U.S.A. 

2 Ostlandring 53, D-25335 Elmshorn, Schleswig-Holstein, Germany 
Version September 11, 2015 

ABSTRACT 

The goal of this investigation is to reinterpret and upgrade the astrometric and other data on comet 
C/1945 XI, the least prominent among the Kreutz system sungrazers discovered from the ground in 
the 20th century. The central issue is to appraise the pros and cons of a possibility that this object is 
— despite its brightness reported at discovery — a dwarf Kreutz sungrazer. We confirm Marsden’s 
(1989) conclusion that C/1945 XI has a common parent with C/1882 R1 and C/1965 Sl, in line with 
the Sekanina-Chodas (2004) scenario of their origin in the framework of the Kreutz system’s evolution. 
We integrate the orbit of C/1882 R1 back to the early 12th century and tlien forward to around 1945 
to determine the nominal direction of the line of apsides and perform a Fourier analysis to get insight 
into effects of the indirect planetary perturbations. To better understand the nature of C/1945 XI, its 
orbital motion, fate, and role in the hierarchy of the Kreutz system, as well as to attempt detecting the 
comet’s possible terminal outburst shortly after perihelion and answer the question in the title of this 
investigation, we closely examined the relevant Boyden Observatory logbooks and identified both the 
photographs with the comet’s known images and nearly 20 additional patrol plates, taken both before 
and after perihelion, on which the comet or traces of its debris will be searclied for, once the process of 
their digitization, currently conducted as part of the Harvard College Observatory’s DASCH Project, 
has been completed and the scanned copies rnade available to the scientific community. 

Subject headings: comets: general - methods: data analysis 


1. INTRODUCTION 

As the most extensive system of genetically related 
comets in existence, the Kreutz sungrazers represent an 
inexhaustible source of research opportunities. Kreutz’s 
(1888, 1891, 1901) celebrated orbital studies described 
the motions of the early bright members, showing that 
they moved about the Sun in similar, extraordinarily 
elongated paths, with orbital periods of up to about 1000 
yr, yet approaching the Sun’s surface to within one solar 
radius at perihelion. Kreutz’s work was followed by many 
more studies in the 20th century, with those by Marsden 
(1967, 1989) standing out as the most important. 

Because all Kreutz sungrazers are fragments of one pro- 
genitor, their orbits’ lines of apsides are nearly perfectly 
aligned; the scatter is only a srnall fraction of 1°, a prod- 
uct of indirect perturbations by the planets, Jupiter in 
particular (Marsden 1967), and of the process of cas- 
cading fragmentation (Sekanina 2002). The research on 
the Kreutz system has recently accelerated explosively 
thanks to vast new evidence provided by imagers on 
board the spacecraft dedicated to the exploration of the 
Sun, especially the coronagraphs C2 and C3 of the Solar 
and Heliospheric Observatory ( SOHO\ see Brueckner et 
al. 1995) and the coronagraplis COR2 and imagers HIl of 
the Solar Terrestrial Relations Observatory's two probes 
( STEREO-A and R; see Howard et al. 2008). Over tlre 
past two decades, these instruments allowed detection, in 
close proximity of the Sun, of thousands of minor Kreutz 
system’s members, referred to hereafter as the dwarf sun- 
grazers, which keep streaming toward the Sun but always 
disintegrate shortly before reaching perihelion. 
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The directions of the lines of apsides of 1600 dwarf 
Kreutz sungrazers derived from their published gravi- 
tational orbits were recently shown by us (Sekanina & 
Kracht 2015; hereafter referred to as Paper 1) to be dis- 
tributed along an arc of 25° (sic!) in the ecliptical lat- 
itude, failing utterly to comply with the condition of 
directional alignment. We determined that this major 
effect was due to a neglected nongravitational accelera- 
tion in the dwarf sungrazers’ motions, which was orders 
of magnitude greater than the nongravitational acceler- 
ations in the motions of the cataloged cornets in nearly- 
parabolic orbits with perihelia a few tenths of AU from 
the Sun or more, topping in exceptional cases tlie Sun’s 
gravitational acceleration. In summary, the preperihe- 
lion disintegration and a very high erosion-driven non- 
gravitational acceleration of tlie orbital motion are two 
fundamental attributes of the dwarf Kreutz sungrazers. 

2. SUNGRAZING COMET C/1945 XI (DU TOIT) 

On 1945 December 11, a cornet was discovered photo- 
graphically by D. du Toit at the Harvard College Obser- 
vatory’s Boyden Station near Bloemfontein, South Africa 
(Paraskevopoulos 1945); nowadays this comet is referred 
to as C/1945 XI. The object moved rapidly toward tlre 
Sun and its brightness at discovery was reported as mag- 
nitude 7. Additional plates were taken at Boyden on the 
following nights until December 15, but the five estimated 
astrometric positions were not communicated until 1946 
January 2, when Cunningham (1946a) used them to com- 
pute three very preliminary parabolic orbits. They indi- 
cated that the comet was apparently a Kreutz sungrazer 
that had passed perihelion five days before the cable was 
sent. Cunningham (1946b) pointed out that his search 
ephemeris, based on one of the three orbits, was uncer- 
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tain “by many degrees.” Published accounts show that 
after Dccember 15 the comet was lost and never seen 
again. Contrary to expectations, it did not becorne a 
brilliant object near and/or after perihelion. Its fail- 
ure to develop a bright headless tail for a limited period 
of time after perihelion — contrary to such sungrazers 
as C/1887 B1 (see the references in Sekanina 2002) and 
C/2011 W3 (e.g., Sekanina & Chodas 2012) — suggests 
that C/1945 XI may have disintegrated already before 
perihclion, as do the dwarf Kreutz sungrazers. This pos- 
sibility raises a question of whether or not tlris object 
was the only dwarf Kreutz sungrazer discovered and re- 
peatedly observed from the ground outside a total solar 
eclipse. 1 To address this issue in detail requires that 
three critical points be answered: 

(1) Can a very liigh nongravitational acceleration be 
conclusively detected in the comet’s orbital rnotion frorn 
the hve Boyden observations? 

(2) Wliy was the comet so bright 17 days before perilre- 
lion? Was it in outburst or was the reported brightness 
grossly overestimated? 

(3) Was the absence of the comet’s relics after perihe- 
lion conditioned on a compelling qualiher or constraint, 
so that the comet’s preperihelion disintegration could be 
subject to doubt? 

We present new evidence in the following sections that 
allows us to connnent on these points and to chart the 
lines of attack in the near future. 

3. THE BOYDEN PHOTOGRAPHIC OBSERVATIONS 

A striking feature of tlie information on the comet’s 
Boyden observations is an extremely slow and protracted 
progress in propagating their results, with an essentially 
complete lack of details. This is most surprising, given 
that C/1945 XI was the hrst Kreutz sungrazer of the 
20th century, after a pause of nearly 60 years. 

An emphatic example of the slow progress is the pub- 
lication of the comet’s astrometry from the Boyden 
plates. After a delayed message of the estimated posi- 
tions (wliose observation tirnes were announced witlr a 
precision to 1 hr!), there was no follow-up report and no 
computation of an improved orbit until 22 years later, 
when Marsden (1967), revealing that the plates were 
measured and reduced by A. G. Mowbray in 1952, de- 
rived several sets of parabolic elements. Strangely, these 
astrometric data were reported by neither Mowbray him- 
self nor L. E. Cunningham, for whorn Mowbray was then 
working (Hockey 2009). In fact, the positions were pub- 
lished only 37 years after tliey were measured and re- 
duced and 44 years after they were taken by the Boyden 
observers (Marsden 1989)! 

The circumstances of tlre comet’s observations at Boy- 
den have never been published. We were especially in- 
terested in the telescopes or cameras employed to make 
these photographic observations, in exposure times used, 
and in the type of tracking (sidereal or on the cornet) of 
the post-discovery photographs. The only information 
learnt frorn the literature was that in the 1940s Boyden 
observers, such as M. J. Bester, discovered their cornets 

1 A dwarf sungrazer C/2008 Ol was detected during a total solar 
eclipse on some exposures after a search based on data from SOHO 
(Pasachoff et al. 2009). C/1882 K1 (Tewfik), discovered ~3.5 hr 
preperihelion at 5.6 Rq (Marsden 1967. 1989) during a total solar 
eclipse on 1882 May 17, may have been a dwarf sungrazer as well. 


witli eitlier the Metcalf 25-crn f/4.9 Triplet refractor or 
the Bache 20-cm f/5.7 Doublet refractor (Cooper 2003, 
2005) while examining plates for irnage quality. 2 

As integral part of the Harvard College Observatory’s 
Digital Accessto a Sky Century@Harvard (DASCH) Proj- 
ect, all Boyden plates are in the process of being dig- 
itized, with the scans gradually made available to the 
scientific community (Simcoe et al. 2006; Grindlay et al. 
2012). 3 The products of this effort will eventually play 
a major role in our quest for information on the comet. 
Of imminent interest to us are the observing logbooks 
on the Harvard website; after some search, we have al- 
ready been able to identify all five plates measured and 
reduced by Mowbray. Three surprises surfaced: (i) the 
comet was not discovered on a plate taken with the 
Metcalf or Baclie telescope; (ii) a never-reported post- 
perihelion search was attempted with two instruments on 
1946 January 8, 11 days after perihelion and the same 
day that Cunningham’s (1946b) ephemeris was issued; 
and (iii) none of the follow-up observations was made by 
the discoverer. 

Table 1 presents the results of our search for the rel- 
evant Boyden plates. For each, the individual columns 
list: the plate number, the UT time of mid-exposure; the 
approximate equatorial coordinates of the plate center 
(converted to the equinox J2000); the local sidereal times 
of the exposure’s start and termination; the resulting ex- 
posure time; the logbook reference; and the observer. 
Three instruments were employed to observe, or searcli 
for, tlie comet: the Metcalf Triplet (plates MF; a plate 
scale of mm _1 ) and two small cameras mounted 

piggyback on the Metcalf, the Cooke 3.8-cm f/8.9 lens 
(plates AM, including the discovery one; a plate scale 
of 610"8 mm _1 ) and the Ross-Fecker 7.5-cm f/7 camera 
(plates RB; a plate scale of 395". 5 mm -1 ). The limit- 
ing magnitudes of the three instruments are listed on 
the Harvard website as, respectively, 17, 13-14, and 15, 
all much fainter than the comet’s reported magnitude at 
discovery. 

All plates employed for C/1945 XI had apparently a 
blue sensitive emulsion (class L; see tlie DASCH website 
footnote) and all three instruments used a plate size of 
20.3 cm (in right ascension) by 25.4 cm (in declination) , 4 
covering fields of 15?1 (Metcaíf), 35?7 (Ross-Fecker), and 
55?2 (Cooke) in angular extent along a diagonal. 

The logbooks show that, of the five observations astro- 
metrically measured and reduced by Mowbray and pub- 
lished by Marsden (1989), tlie ones on December 11 12 
were made with the Cooke lens and only the ones on De- 
cember 13-15 with the Metcalf refractor. The logbooks 
further show that the comet was ~12° from the center 
on the discovery plate, 2°.9 on the December 12 Cooke 
plate, always less than 1°.2 on the Metcalf plates, and 
just about 9°.6 on the December 14 Cooke plate, which 
apparently has never been measured. The question of 

2 These plates were taken on behalf of H. Shapley, Director of the 
Harvard College Observatory, for the purpose of studying southern 
variable stars (Cooper 2003, 2005; van Heerden 2008). 

3 The DASCH Project, at http://dasch.rc.fas.harvard.edu/, 
involves digitization of more than 500 000 Harvard plates and is cur- 
rently in progress; see also http://tdc-www.harvard.edu/plates/. 

4 The website http://tdc-www.harvard.edu/plates/mf/ speci- 
fies that Metcalf plates have generally sizes 8 in by 10 in, 14 in by 
17 in, and, some early ones, 10 in by 12 in. 
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Table 1 

Boyden Plates with Reported Images or Possible Images of Comet C/1945 XI. 


Plate 

number a 

UT time at 
mid-exposure 

Plate center b 

Sidereal time 

Exposure 

time 

(min) 

Logbook 

reference 


R.A. 

Dec. 

start 

stop 

Observer c 

AM25201 d 

1945 Dec. 11.04709 

13 h 36 m 7 

—60°31' 

7 h 26 m 

8 h 56 m 

90 

am45b_0158 

du Toit 

AM 25206 

12.04713 

15 18.7 

-65 22 

7 30 

9 00 

90 

am45b_0158 

Bester 

MF 34987 

13.07210 

16 19.0 

-63 03 

8 40 

9 10 

30 

mf25b_0086 

Britz 

AM 25208® 

14.03336 

16 40.9 

-70 12 

7 43 

8 23 

40 

am45b_0158 

j(?) 

MF 34991 

14.05898 

16 42.9 

-61 12 

8 25 

8 55 

30 

mf25b_0086 

j(?) 

MF 34994 

15.06906 

17 02.6 

-58 39 

8 51 

9 06 

15 

mf25b_0088 

j(?) 

MF 35030 

1946 Jan. 8.06517 

17 17.5 

-47 19 

10 05 

10 50 

45 

mf25b_0094 

Bester 

RB 14184 

8.07729 

17 17.5 

-47 19 

10 25 

11 05 

40 

rbl0_148 

Bester 


a AM plates taken with the Cooke lens; MF plates with the Metcalf telescope; and RB plate with the Ross-Fecker camera. 
b Equinox J2000. 

c The identity of the observer whose abbreviation in the logbooks was J could not be determined; see also Table 12. 
Discovery plate. 

e This plate appears to have never been astrometrically measured and reduced. 


whether or not the comet is located on the plates taken 
on January 8 can only be answered after our orbital anal- 
ysis (see Sections 7.1 and 8.2). 

4. HISTORY OF THE COMET’S ORBIT DETERMINATION 

The sets of orbital elements for cornet C/1945 XI com- 
puted by Cunningham (1946a) were rnuch too uncertain 
to use as a basis of subsequent research. A more accu- 
rate orbit was a gravitational parabolic approximation 
derived by Marsden (1967), who employed the astrome- 
try by Mowbray. Marsden concluded that the December 
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Figure 1. Line of apsides for seven bright Kreutz sungrazers and 
C/1945 XI. The plot of the perihelion longitude against per- 
ihelion latitude B n shows that C/1945 XI, represented by Mars- 
den’s cataloged gravitational orbit, deviates significantly from the 
cluster of the bright sungrazers, whose lines of apsides are closely 
aligned. The asterisk marked Pert. is the apsidal-line position for 
C/1945 XI on the assumption that the comet is a fragment of a 
common parent with C/1882 Rl, the difference being due entirely 
to the indirect planetary perturbations. For six comets the errors 
of the coordinates L^, are smaller than the size of the symbols; 
the estimated errors for C/1880 C1 and C/1945 XI are depicted. 
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14 data point (which we found was measured on a MF 
plate) was inconsistent with the other four positions and 
he omitted it from what he described as the final set 
of orbital elements obtained by least squares. The four 
employed positions were fitted witli a mean residual of 
±2".l, the orbit appeared to be similar to those of C/1882 
R1 and C/1965 Sl, and ever since 1972 it has been listed, 
as the representative orbit for C/1945 XI, in the Cata- 
logue of Cometary Orbits (see Marsden & Williams 2008 
for the most recent edition). 

More recently, this comet’s orbital motion was further 
examined by Marsden (1989). He computed four differ- 
ent gravitational solutions, A D, constraining them to a 
prescribed orientation of the line of apsides and allowing 
the orbit to depart from a parabola. His solutions C and 
D are discussed more extensively in Section 6.1. 5 

In terms of the orientation of the line of apsides, a 
very stable orbital parameter, the representative orbit of 
C/1945 XI is compared in Figure 1 with more reliably de- 
termined orbits of seven bright sungrazers. The errors for 
six of them are smaller than the size of the symbols. For 
C/1945 XI the error is estimated by comparing the scat- 
ter in the apsidal line with tliat in the angular elements, 
combined with the errors of the elements, as published 
by Marsden (1967). The estimated error for C/1880 C1 
is due largely to the uncertainty in the orbital period, as 
investigated by Kreutz (1901). There is a striking dis- 
crepancy between C/1945 XI and the other sungrazers in 
Figure 1, exceeding six standard deviations, thus provid- 
ing grounds for suspecting that the motion of C/1945 XI 
might be — like the motions of the dwarf sungrazers — 
significantly affected by nongravitational forces. 

As of now, it has not been demonstrated conclusively 
whether there exists a purely gravitational orbital so- 
lution for C/1945 XI that simultaneously (i) fits satis- 
factorily at least four of the five available observations; 
(ii) complies with the proper orientation of the line of ap- 
sides; and (iii) is consistent with a plausible osculating 
orbital period. Addressing this issue requires in the first 
place that we know the appropriate values of the quan- 

5 Marsden (1989) remarked that the December 14 observation 
was reconstructed for an arbitrarily chosen observation time from 
the residuals in his previous paper (Marsden 1967), because the 
original data were unfortunately lost. 




4 


Sekanina & Kracht 


tities under (ii) and (iii), which depend on the perturba- 
tions by the planets, Jupiter in particular. We examine 
these perturbation effects next. 

5. INDIRECT PERTURBATIONS BY THE PLANETS 

Because of the nature of their orbits, the Kreutz cornets 
cannot experience a close approach to the planets, in- 
cluding Jupiter. Nevertlicless, the sungrazers’ orbits are 
subjected to indirect planetary perturbations over the 
entire revolution about the Sun and, as a result, show 
limited variations. 

Addressing this problem in some detail, Marsden 
(1967) began his investigation of the effects of indirect 
perturbations on a fragment, separating from its sungraz- 
ing parent at perihelion, by integrating tlre fragment’s 
rnotion over tlre orbital period to the tirne t n of next 
perihelion. First, he considered only Jupiter in a circular 
orbit in the plane of the ecliptic. He applied the equa- 
tions for the variation of arbitrary constants and found 
that in this simplified scenario the fragment’s orbital el- 
ements — the argument of perihelion u>, the longitude of 
the ascending node Í2, the inclination i, and the perihe- 
lion distance q — at time t n depend on Jupiter’s ecliptical 
longitude at this same time, Aj)^), as follows: 

w(Aj) = u> 0 + X u sin(Aj+A 0 ), 

^(Aj) = Ho + Aq sin(Aj+Ao), 

*(Aj) =i 0 + Xi sin(Aj+A 0 ), 

<?(Aj) = qo - X q sin(Aj+A 0 —90°) 

= q 0 +X q cos(Aj+A 0 ), (1) 

where w 0 , ■■■ 1 q 0 are constant terms, X u , ...,X q are 
amplitudes (all taken as positive numbers), and A 0 is a 
constant phase shift. Because of Jupiter’s nonzero orbital 
eccentricity and the deviation of its orbital plane frorn the 
plane of tlre ecliptic, and also because of the indirect per- 
turbations by the other planets, the quasi-periodic vari- 
ations of the sungrazers’ orbital elements are rnore com- 
plex. Marsden indicated that when only Jupiter’s per- 
turbations were accounted for, tlre amplitudes amounted 
to approximately X u = l°.l, Xq, = 1°.4, X, = 0°.3, and 
X q = 0.00039 AU = 0.084 R Q . Marsden furtlrer pointed 
out that his numerical integrations, which took into ac- 
count the perturbations by the planets Jupiter to Nep- 
tune and by Pluto, resulted in expressions similar to 
Equation (1), with the amplitudes equal, respectively, 
to 1°.6, 2°.l, 0°.4, and 0.00046 AU (= 0.099 R Q ). 

Since the longitude and latitude of perihelion, L ïï and 
Bjr, are related to the three angular elements by 

tan(L T —U) =tanw cosí, 

sinBjr =sinwsini, (2) 

the variation of the line of apsides is described, to a first 
approximation, by: 

L n (A 0 ) = L 0 + Xl sin(Aj+A 0 ), 

B n (A 0 ) = B 0 + Xb sin(Aj+A 0 ), (3) 

where 

X L = Xq + (X w cos i 0 — Xi cos w 0 sin B 0 ) sec 2 B 0 , 

Xb = {Xu cot oj 0 + Xi cot t 0 ) tan B 0 . (4) 

Marsden noticed tliat the line of apsides is scarcely af- 


fected by the indirect perturbations. To get a more pro- 
found insight into their infiuence on the orbital rnotion 
of C/1945 XI, we first integrated the orbit of its major 
presumed sibling, C/1882 R1 (see Table 1 of Sekanina & 
Chodas 2004), 6 back in tirne to tlie 1106 perihelion. Next 
we varied the eccentricity incrementally and ran a set of 
orbits, obtained in this faslrion, forward in tirne. Using 
the JPL DE405 ephemeris, we accounted for tlie pertur- 
bations by all eight planets, by Pluto, and by the three 
rnost massive asteroids, as well as for the relativistic ef- 
fect. The increments were adjusted stepwise so that the 
perihelion times ranged between 1939 August 2 and 1951 
August 29 at nearly constant intervals of about 100 days 
each; tlre length of the covered time period slightly ex- 
ceeded Jupiter’s orbital period. Clrecks showed that the 
relative errors accunmlated over an integration period of 
~1600 yr did not exceed ~10~ 10 in the comet’s position 
vector and its velocity vector. 

The integrated perturbations of the six orbital parame- 
ters, ijJ, Í2, i, L^, B n , and q, for tliis period of sungrazers’ 
arrival to perihelion are presented in Figure 2 as a func- 
tion of Jupiter’s longitude Aj at the perihelion times, t n . 
The variations are not periodic in that tlre values of the 
elements at Aj = 0° in 1939 and 1951 differ by 0°.70 in u>, 
by 0°.59 in Q, by 0°.12 in i, by 0.068 R Q in q, by 0°.044 
in Ltt, and by 0°.012 in B n . Frorn Figure 2 it follows 
that the variations for the last four parameters are rea- 
sonably well approximated by the thick sine curve, while 
those for u> and íí are seen to deviate quite significantly, 
especially in the range of Jupiter’s longitudes frorn 30° to 
150°. For at least tliese two elements the scenario based 
on Jupiter in a circular orbit in the plane of the ecliptic 
is clearly inadequate. 

The perturbation variations in the perihelion longitude 
Ltt are a problem. Equation (3) suggests that the curve 
should be in phase with the variations in the argument 
of perihelion, the longitude of the ascending node, and 
the inclination, but Figure 2 demonstrates that it is not. 
Closer inspection shows that the reason for this discrep- 
ancy is the fact that the predicted amplitude given by 
Equation (3) is very close to zero (or, in fact, slightly 
negative) and that the actual variations in L n in Figure 2 
are determined by the second order terrns tliat have been 
neglected in Equation (3). For the perihelion latitude B n 
this is not the case and its variations are indeed in phase 
with u>, Í2, and i. 

The numerical integration shows tlrat the net ampli- 
tudes of the orbital parameters amount to 1°.52 in lo, 
1°.86 in íl, 0°.39 in i, 0°.ll in L n (which at a latitude of 
~35° represents an arc of 0°.09), 0°.07 in B n , and 0.00051 
AU or 0.110 R Q in q. Tlrese are typically within about 
10% of the amplitudes given by Marsden for the four el- 
ements from Equation (1) above. The overall amplitude 
for the line of apsides is 0°.ll, so that the range, 0°.22, is 
generally consistent with the apsidal difference of 0°.30 
between C/1882 R1 and C/1843 D1 accumulated over a 
span of two orbital periods (or about 15-17 centuries), 
suggested for tlre age of the Kreutz systern by Sekanina 
& Chodas (2004, 2007). 

6 This presumption is based both on Marsden’s (1967, 1989) and 
Sekanina & Chodas’ (2004) conclusions and on the results of our 
own experimentation with the orbit of C/1945 XI, which showed 
that it was much easier to align its apsidal orientation with that of 
C/1882 R1 than C/1843 Dl, the other potential major sibling. 
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Figure 2. Plots, against the ecliptical longitude of Jupiter Aj at the time of a respective sungrazer’s perihelion, of the argument of 
perihelion lj, the longitude of the ascending node the inclination i, the perihelion longitude and latitude B^ (all equinox J2000), and 
the perihelion distance q of the siblings of comet C/1882 R1 separated at the 1106 perihelion, which returned to perihelion between 1939 
August 2 and 1951 May 18. The sungrazer C/1945 XI was at Aj = 195°.23. The thick sine curves show the simple approximations of the 
type expressed by Equation (1), the thinner curves are the much better approximations by the Af-harmonic Fourier polynomials (N = 4 for 
all but the last curve, for which N = 2; see Table 3). Note that the scales for and are four times wider than those for u ;, íí, and i. 
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Table 2 

Cataloged Orbit of C/1945 XI and Orbit of the Hybrid 
(Equinox J2000). 


Orbital element 

C/1945 XI 
(cataloged a ) 

Hybrid b 

(perturbed) 

Perihclion time t n (1945 TT) 

Dec. 27.9652 

Dec. 27.9801 

Argument of perihelion cj 

72°,0619 

67°.8518 

Longitude of ascending node O 

351°.2006 

345°.3759 

Orbit inclination i 

141°,8734 

141°,5608 

Perihclion distanc.e q { ^ U \ 

0.007516 

0.007123 

l (u©) 

1.6147 

1.5302 

Orbit eccentricity 

1.0 

0.99992634 

Semimajor axis (AU) 

oo 

96.7 

Orbital period (yr) 

oo 

951 

Longitude of perihelion 

283°.57 

282°.84 

Latitude of perihelion 

+35°.97 

+35°. 16 

Osculation epoch (1945 TT) 

(none) 

Dec. 26.0 


a See Marsden 8z Williams (2008). 

b Orbit of C/1882 R1 integrated to perihelion time of C/1945 XI. 


We now constrain the orbital properties of C/1945 XI 
by its arrival time at perihelion, 1945 December 28.0 TT. 
The predicted orientation of the comet’s line of apsides 
is determined by the ecliptical coordinates L n = 282°. 84 
and B n = +35°. 16 (Equinox J2000), deviating frorn the 
orientation of the apsidal line of C/1882 R1 by 0°.ll. 
The longitude of Jupiter at the tirne of the comet’s per- 
ihelion is Aj = 195°.23. In addition, the computations 
show that the comet’s osculating semimajor axis should 
have been 96.7 AU and the orbital period 951 yr, longer 
than the actual tirne span since 1106. This orbital period 
is in excellent agreement with the value tlrat Sekanina & 
Chodas (2004) obtained for this cornet in their liierarchy 
model of the Kreutz system (see their Table 12). 

To illustrate the influence of the indirect perturbations 
specifically on the orbit of C/1945 XI, we compare Mars- 
den’s (1967) cataloged parabolic orbit (cf. Marsden & 
Williams’ 2008) with tlie orbit of what we call a hybrid 
one that C/1882 R1 would have had, if it arrived at 
perihelion at the tirne C/1945 XI did. The two orbits 
differ by several degrees in the argument of perihelion 
and the longitude of the ascending node, and by almost 
0.1 Rq in the perihelion distance. 

Returning now to Figure 2, one could tliink of the dif- 
ferences between the sequence of the points, computed 
by numerical integration of the orbits, and the thick sine 
curves that are intended to ht them, as deformations and, 
accordingly, employ the superposition principle to miti- 
gate the discrepancies. The effects of the perturbations, 
even though not periodic on a scale of Jupiter’s orbital 
period, can nonetheless be expressed as a combination 
of periodic variations in terms of an Wharmonic Fourier 
polynomial. Figure 2 shows the results of this fitting, 
witli N not exceeding 4, as the lighter curves; obviously, 
the improvement over a sirnple sine curve is considerable. 
Thus, calling 3? any of the six parameters in Equations 
(1) and (3), we can write it in the forrn: 

N 

5ft(Aj) = % + E [oae,fc sin(fcAj) + b^ k cos(fcAj)] 

fc=i 

N 

= Sfto+E^sin(Mj+A K , fc _i), (5) 

fc=i 
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HARMONIC N OF FOURIER POLYNOMIAL 

Figure 3. Mean residuals <j of the Fourier polynomials, fitted to 
the variations in the longitude of the ascending node, the argument 
of perihelion, and the inclination during Jupiter’s orbital period 
1939-1951, as a function of the polynomials’ harmonic N. The 
minimum residual for all three orbital elements is reached at N = 4. 

where 

Vjf,fc = \J a%fc k + b^ k , 

A= arctan Í—L\ . (6) 

V a 5R.fc/ 

It is apparent that for N = 1, Equation (5) emulates 
Equations (1) and (3), with Asr,o being always equal 
to Ao for !ft = (w,íl,*, L n , B n ) and Aj?,o = Ao + 90° for 
5ft = q, and with Xj},! = X/r for any orbital parameter. 

The coefhcients asR,fc and bsR,k are computed by least 
squares from the values of 5ft(Aj) derived from the or- 
bit integration runs as a function of Jupiter’s longitude 
Aj (t n ). The appropriate number of terms N = N m - ln > 1 
is determined by a minimum mean residual. As the 
number of Fourier coefficients in the equations of con- 
dition increases with N, the number of the degrees of 
freedom drops, thus causing the mean residual a to in- 
crease. A byproduct is an increasing uncertainty of 
the Fourier coefficients in tlie (meaningless) polynomi- 
als with N > N min . As examples, the plots of the mean 
residuals for lo, O, and i vs N are displayed in Figure 3. 

The Fourier solutions for the six orbital parameters in 
Figure 2 are listed in Table 3. For eacli parameter the 
first line provides the results of tlie approximate solution, 
fit A, by applying Eciuation (1) or (3). The second and 
following lines refer to the fitted Fourier polynomial, F 'n ■ 
The mean residual a is always followed by the constant 
term 3?o and by information on each harmonic k in the 
case of the V-harmonic Fourier polynomial: the harmonic 
amplitude XsR k and phase Ajf,fc_i. 

We note that the amplitudes of all harmonics are de- 
termined with an error not exceeding ±35%, but mostly 
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Table 3 

Fourier Polynomials Fitting the Dependence of Indirect Planetary Perturbations of Orbital Elements of 
Kreutz Sungrazers in 1939-1951 on Ecliptical Longitude of Jupiter. 


Orbital element 

Fit a 

Mean 

residual 

Constant 

term b 

Harmonic 
k = l,... ,N 

Harmonic 

amplitude 

Harmonic 

phase 

Argument of perihelion, lj 

A 

±0°269 

688454 ± 08057 


18397 ± 08083 

35682 ± 383 


f 4 

±0.084 

68.455 ± 0.018 

i 

1.386 ± 0.026 

355.9 ± 1.0 





2 

0.272 ± 0.026 

202.1 ± 5.3 





3 

0.183 ± 0.025 

57.3 ± 8.0 





4 

0.105 ± 0.025 

319 ± 14 

Longitude of ascending node, 

A 

±0.338 

346.214 ± 0.072 


1.731 ± 0.104 

359.1 ± 3.3 


f 4 

±0.113 

346.216 ± 0.024 

1 

1.718 ± 0.035 

358.8 ± 1.1 





2 

0.338 ± 0.035 

205.8 ± 5.8 





3 

0.229 ± 0.034 

59.0 ± 8.7 





4 

0.134 ± 0.035 

323 ± 15 

Orbit inclination, i 

A 

±0.061 

141.709 ± 0.013 


0.364 ± 0.019 

357.0 ± 2.8 


f 4 

±0.023 

141.710 ± 0.005 

1 

0.361 ± 0.007 

356.7 ± 1.1 





2 

0.059 ± 0.007 

198.6 ± 6.7 





3 

0.042 ± 0.007 

52.0 ± 9.6 





4 

0.024 ± 0.007 

320 ± 17 

Longitude of perihelion, 

A 

±0.0202 

282.9179 ± 0.0043 


0.0968 ± 0.0060 

104.6 ± 3.9 


f 4 

±0.0087 

282.9183 ± 0.0019 

1 

0.0906 ± 0.0026 

104.1 ± 1.7 





2 

0.0222 ± 0.0026 

299.2 ± 6.9 





3 

0.0077 ± 0.0027 

143 ± 20 





4 

0.0082 ± 0.0027 

29 ± 18 

Latitude of perihelion, 

A 

±0.0169 

±35.1860 ± 0.0036 


0.0616 ± 0.0052 

354.4 ± 4.6 


f 4 

±0.0030 

±35.1859 ± 0.0006 

1 

0.0608 ± 0.0009 

353.8 ± 0.9 





2 

0.0185 ± 0.0009 

194.6 ± 2.7 





3 

0.0111 ± 0.0009 

66.3 ± 4.7 





4 

0.0055 ± 0.0009 

321.4 ± 9.3 

Perihelion distance, q (Rq) 

A 

±0.0156 

1.6200 ± 0.0034 


0.0987 ± 0.0048 

84.0 ± 2.9 


f 2 

±0.0132 

1.6183 ± 0.0028 

1 

0.0905 ± 0.0039 

83.2 ± 2.6 





2 

0.0173 ± 0.0040 

325 ± 13 


a A = approximation by Equation (1) or (3), based on the assumption of a single planet in circular orbit in the plane of the ecliptic; 
Fív = A/-harmonic Fourier polynomial that provides the best fit (minimum mean residual). 
b Equinox J2000. 


less than ±25%. With one borderline exception the 
amplitudes decrease with increasing harmonic, as ex- 
pected. The amplitude of the first harmonic is always 
only sliglrtly smaller than the net amplitude derived di- 
rectly from tlie results of numerical integration of the or- 
bit. The phase angles of all harmonics for the argument 
of perihelion, the longitude of the ascending node, the in- 
clination, and the latitude of perihelion are remarkably 
close to one another. The first harmonic for the perilre- 
lion distance is shifted, as expected, by alrnost exactly 
90° relative to those of w, ÍI, i, and B^. Only the lon- 
gitude of perihelion L n appears to be out of phase with 
any of the other orbital parameters. 

The least-squares procedure also allows us to generalize 
Equation (5) thus: 

N 

5i(Aj) = 3? 0 ± ^2 -AsR.fe sin(cfcAj±Aïf ifc „i), (7) 

fc=í 

where c^’s are any positive numbers, not restricted to 
integers. This generalization offers us an opportunity to 
relate the coefficients Ck (k > 1) to the orbital periods of 
other planets, Pk, so that Ck = Pj/Pk- Without provid- 
ing any details, we point out that our search for least- 
squares solutions of this type proved rnuch less successful 
in fitting the quasi-periodic variations in the sungrazers’ 
orbital elements tlian did the Fourier polynomials. 


6. ORBITAL COMPUTATIONS BASED ON THE BOYDEN 
OBSERVATIONS FROM DECEMBER 11 15 

If comet C/1945 XI is closely related to C/1882 Rl, 
the hybrid’s orbit in Table 3 suggests tliat the line of 
apsides be described by L n = 282?84 and B n = ±35?16 
(Equinox J2000), as plotted in Figure 1 (marked Pert.), 
and tliat the osculating semimajor axis at perihelion be 
equal to a = 96.7 AU. Table 3 also shows that the oscu- 
lating elements sliould differ from those of C/1882 R1 by 
~2° in u! and f2, by almost 0°.5 in i, and by more than 
0.1 Rq in q. Below we investigate a range of orbital solu- 
tions, gravitational and nongravitational, offered by the 
observations made at Boyden on December 11 15. 

As measured and reduced by Mowbray and listed for 
the equinox B1950 by Marsden (1989), the five astro- 
metric positions are for the equinox J2000 presented in 
Table 4. Because on the first two nights the comet was 
imaged with tlie Cooke lens, while on the three remain- 
ing nights with the Metcalf refractor, we assign the posi- 
tional data weights that are inversely proportional to the 
plate scales of the two instruments. Comparison of the 
observation times in column 2 of Table 4 with tliose in 
colunm 2 of Table 1 shows that the first ones are system- 
atically too early by 0.00021 to 0.00022 day, or 18—19 s. 
The nature of the difference is unknown; it cannot come 
from a longitude discrepancy, because it is equivalent to 
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Table 4 

Measured and Reduced Boyden Observations of Comet C/1945 XI (Equinox J2000). 


Ref. 

No. 

Observation time 
1945 (UT) 

R.A. 

Dec.. 

Assigned 

weight 

Instrument 

employed 

Plate 

number 

From plate center 

distance 

RA. 

1 

December 11.04687 

15 h 13 m 24?42 

—65°23 43"4 

1.0 

Cooke 

AM 25201 

11394 

12434 

2 

12.04691 

15 43 21.25 

-64 09 13.1 

1.0 

Cooke 

AM 25206 

2.89 

68.0 

3 

13.07189 

16 11 04.21 

-62 29 32.5 

3.6 

Metcalf 

MF 34987 

1.07 

300.5 

4 

14.07000 a 

16 34 56.77 

-60 33 16.5 

3.6 

Metcalf 

MF 34991 

1.16 

302.6 

5 

15.06885 

16 55 41.67 

-58 22 47.2 

3.6 

Metcalf 

MF 34994 

0.95 

285.7 


a Position reconstructed by Marsden for an arbitrarily chosen time on this date (see footnote in Section 4). 


a distance of as much as ~500 meters. The UT times in 
Table 1, determined from the sidereal tirnes of the be- 
ginning and end of the exposures, are in fact UTl times 
and have independently been checked; they never left an 
unexplained difference of more than ~1 s. Nevertlieless, 
for the sake of comparison with Marsden’s (1967, 1989) 
results, we retain for this exercise his UT tirnes. 

6.1. Gravitational Solutions 

Since the observed orbital arc is too short to reliably 
compute the eccentricity, it needs to be determined frorn 
other considerations. As already mentioned above, the 
constraint is here provided by adopting the semimajor 
axis of tlre hybrid’s osculating orbit, 96.7 AU at the epoch 
of 1945 December 26.0 TT, for all solutions that follow. 
We begin with a gravitational orbit to fit all five astro- 
metric positions in Table 4. 

Deriving orbits B and C and comparing tlrem with or- 
bit A, Marsden (1989) demonstrated that elliptical solu- 
tions fitted tlie five observations of C/1945 XI rnuch bet- 
ter than a parabolic approximation. Yet, neither orbit, 
B or C, provides a satisfactory fit. One possible contrib- 
utor to this problem could be the constraint introduced 
to satisfy tlre prescribed direction of the line of apsides 
(for which Marsden took L v = 282°. 7, B v = +35°. 2 af- 
ter conversion to Equinox J2000), the sarne issue that 
we encountered with the dwarf sungrazers of the Kreutz 
system (Paper 1; see also Section 1). 

Comparison of orbit C, the rnore realistic of tlie two 
sets in terms of the comet’s orbital period, with our opti- 
mum gravitational fit to the five observations is presented 
in Table 5. Marsden’s (1989) residuals frorn orbit C in 
this table have a tendency to become more positive to- 
ward both ends of the observed arc in right ascension and 
to get more negative with tirne in declination. We tried 
to emulate the residuals frorn orbit C, employing the set 
of truncated elements in Marsden’s (1989) Table V, but 
succeeded to achieve this within about 4" only in decli- 
nation. Although we allowed for the rounding off of the 
perihelion time, the residuals in right ascension carne out 
systematically more negative by about 12" and all except 
the last one were negative. This effect is due apparently 
to the elements’ truncation. 

The residuals from our Best Fit solution, although 
clearly better than those from orbit C, are not quite 
satisfactory either, suggesting tliat the first observation 
may be inferior. That should not be surprising, consid- 
ering that the comet was rather far from the plate center 
on this discovery exposure (Section 3). Unfortunately, 
the line of apsides deviates from the proper direction by 


0°.571, which is unacceptably large. Although tlie in- 
dividual angular elements of the Best Fit solution are 
burdened with errors of up to several tenths of a degree, 
the correlation among them draws the uncertainty in the 
orientation of the line of apsides down to only hundredths 
of a degree. 

Next, we compare Marsden’s (1989) orbit D witli the 
hybrid’s orbit frorn Table 2. As with orbit C, we first 
tested whether the low-precision orbit D from Marsden’s 
Table V is able to reproduce the residuals in his Table VI. 
Since this exercise did not rneet with success, we decided 
to reconstruct the high-precision version of orbit D, not 
published by Marsden, by subtracting the residuals in 
his Table VI from the observed astrometric positions, re- 
taining his constraint on the reciprocal semimajor axis 
(0.01016 AU -1 ). This approacli was successful; the re- 
sulting orbit D is sliown in Table 6 to be exceedingly 
similar to tlie hybrid’s orbit, the two sets of elements dif- 
fering mostly in the fourth or higher significant digit. In 
particular, we note that the lines of apsides agree to bet- 

Table 5 

Marsden’s (1989) Orbit C and Our Optimum Gravitational 
Solution to Fit All Five Boyden Observations of 
C/1945 XI (Equinox J2000). 


Orbital element 

Orbit C 

Best Fit 

Perihelion time t n (1945 TT) 

Dec. 27.976 

Dec. 27.982 

Argument of perihelion lj 

68°. 10 

70°.99 

Longitude of ascending node Lí, 

345°.54 

349°.50 

Orbit inclination i 

141°.60 

141°.91 

Perihelion distance q { / 

0.00699 

1.502 

0.007244 

1.556 

Orbit eccentricity 

0.99993010 

0.99992506 

Semimajor axis (AU) 

100 

96.67 

Orbital period (yr) 

1000 

950.5 

Longitude of perihelion 

282°. 70 

283°. 14 

Latitude of perihelion 

+35°.20 

+35°.68 

Osculation epoch 1945 (TT) 

(none) 

Dec. 11.0 

Mean residual a 

±11395 

±6317 


Distribution of Residuals a , O — C 


Time of 
observation 

1945 (UT) 

Orbit C 

Best Fit 

R.A. 

Dec. 

R.A. 

Dec.. 

Dec. 11.04687 

+536 

+2234 

+1332 

+830 

12.04691 

+6.1 

+5.3 

+6.1 

+2.5 

13.07189 

-14.8 

-11.8 

-6.6 

-5.6 

14.07000 

+ 10.3 

-1.9 

+7.8 

+8.0 

15.06885 

+ 15.6 

-10.9 

-2.7 

-3.2 

a From unweighted observations for orbit C, weighted for Best Fit. 
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Table 6 

Marsden’s (1989) Orbit D and Hybrid’s Orbit (Equinox J2000). 


Orbital element 

Orbit D 

Hybrid’s orbit 

Perihelion time t n (1945 TT) 

Dec. 27.9798 a 

Dec. 27.9801 

Argument of perihelion uj 

67°. 8356 

67°.8514 

Longitude of ascending node Ll 

345°.3571 

345°.3753 

Orbit inclination i 

141°,5553 

141°,5607 

Perihelion distance q l / 

1 (%) 

0.007129 

1.5315 

0.007123 

1.5302 

Orbit eccentricity 

0.99992756 

0.99992632 

Semimajor axis (AU) 

98.42 

96.67 

Orbital period (yr) 

976.4 

950.5 

Longitude of perihelion b , 

282°.84 

282°.84 

Latitude of perihelion b , B n 

+35°.16 

±35°.16 

Osculation epoch 1945 (TT) 

Dec. 11.0 

Dec. 11.0 

Mean residual c < 

f ±7'(63(M) 

( ±8"17(SK) 

±13'Í75 


Distribution of Residuals c , O — C 


Time of 
observation 
1945 (UT) 

Orbit D (M) 

Orbit D (SK) 

Hybrid’i 

s orbit 

R.A. 

Dec.. 

R.A. 

Dec. 

R.A. 

Dec. 

Dec. 11.04687 

±10'(8 

+5"1 

+12'(3 

+6'(6 

+26'(4 

—6'Í3 

12.04691 

-2.5 

-2.8 

-1.1 

-1.1 

+13.7 

-12.1 

13.07189 

-15.2 

-11.4 

-14.0 

-9.6 

+1.2 ■ 

-18.7 

14.07000 

+4.3 

+5.6 

+5.3 

+7.6 

+20.4 

0.0 

15.06885 

+2.8 

+2.5 

+3.6 

+4.5 

+18.4 

-1.8 


a To fit weighted observations, t n — Dec. 27.9800. 
b These ecliptical oordinates define the reference line of apsides. 
c Entry (M) and columns 2—3 from unweighted observations; entry 
(SK) and columns 4—5 from weighted observations; residuals from 
hybrid’s orbit are from weighted observations. 

ter than 0°01. Tlris correspondence is not fortuitous, as 
both sets of elements represent perturbed versions of the 
orbits of two comets observed 83 years apart, wliose rno- 
tions in space were virtually identical, as already pointed 
out by Marsden (1967): C/1882 Rl, used by us to derive 
the hybrid’s orbit; and C/1965 Sl, employed by Marsden 
(1989) to derive orbit D. 

Table 6 also lists three sets of residuals: those left by 
orbit D frorn tlie unweighted observations, in columns 2 
and 3 [taken frorn Marsden’s (1989) Table VI and rnarked 
M]; tlrose left by orbit D from the weighted observations, 
as derived by us' (rnarked SK), in columns 4 and 5; and 
those left by the liybrid’s orbit from the weighted obser- 
vations, in columns 6 and 7. 

Even though the residuals in Table 6 are fairly large, 
the ones left by orbit D are rnuch better than the strongly 
systematic ones left by the hybrid’s orbit. Frorn tliis com- 
parison as well as from the difference in the mean residual 
we infer that the orbital evolution of cornet C/1945 XI, 
as a fragment of its connnon parent with C/1882 R1 and 
C/1965 Sl, was apparently more sinhlar to the orbital 
evolution of the latter than the former. This is an im- 
portant though tentative conclusion, consistent with the 
fragmentation hierarchy of the Kreutz systern proposed 
by Sekanina & Chodas (2004). According to their evolu- 
tionary rnodel, C/1882 Rl, C/1965 Sl, and tlre precur- 
sor to C/1945 XI separated from their comrnon parent, 
possibly X/1106 Cl, at tlie sarne time, some 18 days 

' This fit required that the perihelion time be increased by 0.0002 
day relative to the tabulated value; otherwise the residuals would 
be slightly greater. 
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after perihelion. Relative to C/1882 Rl, the precursor 
moved in the same direction as, and with a velocity only 
about 20% lower than, C/1965 S1 (see also Sekanina & 
Chodas 2002). The subsequent separation of C/1945 XI 
from its precursor around 1700 AD notwithstanding, the 
comet’s orbit in 1945 should indeed resemble the orbit of 
C/1965 S1 to a greater degree than that of C/1882 Rl. 

Since orbit D and the hybrid’s orbit are so very similar, 
yet the residuals they leave substantially differ, we tested 
whether this effect is due to the minor differences in the 
angular elements or the orbital dimensions. We replaced 
the semimajor axis 98.4 AU with 96.7 AU and noted that 
the residuals frorn orbit D in Table 6 did not change; the 
new orbit, D', is now: t^ = 1945 December 27.9803 TT, 
q = 0.007126 AU = 1.5309 R & e = 0.99992630, and for 
the equinox of J2000, ui = 67°.8377, U = 345°.3558, and 
i = 141°.5570. The residuals in Table 6 reflect thus entire- 
ly the differences of up to 0?02 in the angular elements. 

Before turning to nongravitational solutions, we take 
notice of a possibility tlrat one of the five astrometric 
observations is inferior and sliould be discarded before 
computing any orbit. Accordingly, we are now going to 
search for solutions that could fit four observations. Such 
solutions will be referred to as where Q stands for 

gravitational and {3} is a progression of reference num- 
bers of the observations that such solutions are based on; 
these numbers are listed in column 1 of Table 4. The Best 
Fit solution presented in Table 5 can now be referred to 
as (/{1,2, 3,4, 5}. As already mentioned above, the dis- 
covery position may be inferior, so that one of tlre tested 
four-observation solutions is Q{ 2,3,4, 5}. The residuals 
from orbit D in Table 6 suggest that the third observa- 
tion may be even worse than the hrst, raising interest in 
the solution Q{ 1,2,4, 5}. Finally, the derivation of the 
cataloged parabolic orbit implied that the fourtlr obser- 
vation failed to fit that solution (Marsden 1967), hence a 
need to test Q{ 1, 2,3, 5}. The results of all three of tliese 
solutions are presented in Table 7. No orbits are being 
computed based on three observations only, as tlrese are 
regarded for our purposes as meaningless. 

Since a semimajor axis a between 96.7 AU and 98.4 AU 
was shown to rnake no difference, tlre four-observation 
solutions below are subjected to the same constraint as 
the Best Fit solution in Table 5, narnely, a = 96.7 AU 
at an osculation epoclr of 1945 December 26.0 TT. The 
solutions listed in Table 7 differ from one another consid- 
erably. The orbits Q{ 2,3,4, 5} and Q{ 1, 2,3, 5} are both 
unacceptable on account of their large offsets from the 
reference apsidal direction (given by the liybrid’s orbit 
and Marsden’s orbit D). The set Q{ 2,3,4, 5} is in addi- 
tion handicapped by the very large residuals left by the 
December 12 observation. Even cursory inspection shows 
that the orbit Q{ 1, 2,4, 5} is by far the most promising of 
the three, implying that it is the December 13 observa- 
tion that is inferior. The offset from the reference line of 
apsides, slightly less than 0°.l, is reasonably low though 
not completely satisfactory. It is somewhat surprising 
that the apparently inferior astrometric position is one of 
those taken with the Metcalf refractor. We woulcl expect 
that the discarded data point sliould be one of the first 
two observations, rnade with the Cooke camera. Under 
the circumstances, the solution Q{ 1,2,4, 5} can at best 
be regarded as marginally acceptable to approximate the 
orbital motion of C/1945 XI. 
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Table 7 


Comparison of Three Gravitational Solutions Based on 
Four Observations of C/1945 XI (Equinox J2000). 


Orbital element 


G{2, 3,4, 5} 

0(1, 2, 4, 5} 

0(1, 2, 3, 5} 

Perihelion time Ítt 

(1945 TT) 

Dec. 27.986 

Dec. 27.979 

Dec. 27.983 

Argument of perihelion uj 

65°.98 

67°. 43 

-4 

. o 

b 

o 

Longitude of ascending node Ll 

342°.84 

344°. 82 

353°.52 

Orbit inclination i 


141°.37 

141°.50 

142°,07 


0.007004 

0.007120 

0.007415 



1.505 

1.530 

1.593 

Orbit eccentricity 


0.99992755 

0.99992635 

0.99992329 

Longitude of perihelion 

282°.54 

282°.80 

283°. 49 

Latitude of perihelion 

+34°. 77 

+35°.09 

+36°.22 

Reference apsidal line’s offset a 

0°.462 

0°.078 

1°.184 

Osculation epoch 1945 (TT) 

Dec. 11.0 

Dec. 11.0 

Dec. 11.0 

Mean residual (weighted) 

±5"11 

±1"79 

±1"06 


Distribution of Residuals b , 

O-C 


Time of 

G{'2, 3, 4, 5} 

0(1,2,4,51 

0(1, 2, 3, 5} 






1945 (UT) 

R.A. Dec. 

R.A. 

Dec. 

R.A. Dec. 

Dec. 11.04687 


+5"3 

+1"1 

—3"3 +3"3 

12.04691 

+25'Í2 +12"0 

-8.2 

-7.4 

+4.9 +2.7 

13.07189 

-5.1 -4.3 


. +0.1 -0.8 

14.07000 

+4.3 +5.5 

+0.1 

+ 1.3 


15.06885 

1 

1 

co 

+0.1 

-0.8 

-0.2 +0.3 


a Reference line of apsides is defined by the hybrid’s and H-. 
^From weighted observations. 


In any case, it is wortli examining whether any mean- 
ingful refinement of the orbit, an apsidal line in particu- 
lar, can be achieved by incorporating a nongravitational 
acceleration, which, if comparable to those affecting the 
motions of tlre dwarf sungrazers, rnight be detectable 
over a span of four days. Besides, suclr alternative so- 
lutions should prove beneficial to estimating a range of 
uncertainties in the comet’s orbital rnotion (Section 7). 

6.2. Standard Nongravitational Solutions 

In order to find out whether there is evidence for non- 
gravitational effects in the motion of comet C/1945 XI 
in the rneager set of astrometric data available, we be- 
gin with a standard formalism of Marsden et al. (1973), 
based on a water-ice sublimation model. The erosion- 
driven nongravitational accelerations in the three cardi- 
nal directions, i.e., in the radial (away frorn the Sun), 
R, transverse, T, and norrnal, N, directions of a right- 
handed RTN coordinate systern that is referred to the 
comet’s orbital plane, are in this formalism expressed by 

a R (r) Ax 

a T (r) = A 2 -gstd(r), (8) 

_a N (r)J [^ 3 . 

where g s td(r) is the standard law employed by Marsden et 
al. (1973), approximating the sublimation rate of water 
ice frorn an isothermal spherical nucleus and normalized 
so that g s td(l AU) = 1, 

( \ —m r / \n- —k 

1+ (^~) ' ( 9 ) 

Here ?u = 2.15, n = 5.093, nk = 23.5, a scaling distance 
ro = 2.808 AU, and a normalization constant a = 0.1113. 


The parameters A\, A 2 , are, respectively, the radial, 
transverse, and normal components of the nongravita- 
tional acceleration at 1 AU from the Sun; their units 
usually employed in orbital studies are 10 -8 AU day -2 . 
Since tlre vector of the acceleration due to tlie generally 
sunward-directed sublimation points away frorn the Sun, 
physically meaningful values of Ai should always be pos- 
itive. The law (9) lias over the past 40 years served ad- 
mirably in countless orbit-determination applications to 
comets with perihelia typically several tenths of AU. 

The numerical procedure that we apply next follows 
the rnethod described in Paper 1. Briefly, since the di- 
rection of tlie line of apsides was found to be a function of 
the nongravitational acceleration, a minimum offset from 
the reference, or nominal, apsidal line (defined here by 
Lj r and B n for the hybrid’s orbit or orbit D; see Table 6) 
serves as the constraint tliat determines the nongravita- 


tional acceleration’s most probable magnitude. 

Similarly to our notation for the gravitational runs, 
we refer to tliese nongravitational solutions by A/’ S ^(A}, 
where {3} is again a progression of reference numbers of 
the employed observations from Table 4, while de- 
notes, on the one hand, Marsden et al.’s (1973) formalism 
of accounting for the nongravitational effect, witli the 
standard law (9) describing the variations with heliocen- 
tric distance; and, on the other hand, one of the two ver- 


sions applied: either solving for the radial component of 
the acceleration with the parameter A\ when X = R; or 


for the norrnal component with A% when X = N. Wlien- 


ever we attempted to solve for both components, the run 
aborted. We recall tliat the radial component (A\ > 0) 
dominates the magnitude of the nongravitational accel- 
eration in the motions of the cataloged cornets in nearly- 
parabolic orbits (Marsden & Williams 2008), while the 









Comet C/1945 XI (du Toit) 

Table 8 

Comparison of the Nongravitational Parameters 
and Test Results for Eight Orbital Solutions 
Based on Marsden et al.’s formalism. 


N ongr avit at ional 
solution 

Parameter a 

A\ or As 

Apsidal- 

line b 

offset 

Mean 

residual 

(weighted) 

W S «{1,2,3,4,5} 

-0.82 4=0.17 

07116 

±6?42 

W«{2,3,4,5} 

+0.67 4=0.18 

0.118 

±5.05 

^{i.2,4,5} 

+0.113 + 0.019 

0.013 

±1.73 

JV«{1,2,3,5} 

-1.70 + 0.40 

0.270 

±3.03 

^11,2,3,4,5} 

-1.87 + 0.23 

0.073 

±6.68 

A/td* (2, 3, 4, 5} 

+1.54 + 0.30 

0.085 

±4.88 

^{1,2,4, 5} 

+0.234 + 0.018 

0.006 

±1.72 

^{1,2,3, 5} 

-4.18 + 0.65 

0.212 

±2.61 


a In units of 10 -5 AUday —2 . 

b Offset from direction of the reference line of apsides (Table 6). 


normal component (A 3 ) was found to contribute signifi- 
cantly to the magnitude of the acceleration that affects 
the motions of the dwarf Kreutz sungrazers (Paper 1). 

Table 8 summarizes the most important results frorn 
the runs based on Marsden et al.’s (1973) standard for- 
malisrn. Altogether we computed eight sucli solutions for 
four different observational sets {3}, the sarne ones as be- 
fore: . ^{1,2,3,4,5}, ^{2,3,4, 5}, AÍ™{1, 2,4, 5}, 
.A/^dU 2,3, 5}; and similarly with N instead of R. Both 
the R and the N versions of the A(,^{1, 2, 3,4, 5} and 
^«{1,2,3, 5} runs are inferior to their respective grav- 
itational solutions (Tables 5 ancl 7) in terrns of the rnean 
residual: ±6"42 (R) and ± 6"68 (N) against ±6"17 in 
the case of the {1, 2,3,4, 5} solution; and ±3".03 (R) and 
±2".61 (N) against ±1".06 in the case of tlre {1,2,3,5} 
solution. In addition, for these two sets of solutions the 
R -version parameters A\ come out to be negative and 
therefore physically meaningless. 

Even thougli the set {2,3,4, 5} avoids the misfortunes 
of the sets {1,2, 3,4, 5} and {1,2,3, 5}, its offsets frorn 
the reference line of apsides fail to improve over the 
corresponding offsets from {1,2,3,4, 5}; in addition, the 
residuals left by tlre December 12 observation are com- 
parable to tliose frorn the equivalent gravitational solu- 
tion (Table 7) and entirely unacceptable. In summary, 
the standard nongravitational solutions based on the set 
{1, 2,4, 5} are by far the best, in terrns of botlr the offset 
frorn the reference apsidal line and the rnean residual. 
Also, both A\ and A3 of the {1,2,4, 5} solutions come 
out in Table 8 to be almost or just about one order of 
magnitude smaller than those of the poor solutions. 

As the final connnents we note that (i) relatively to 
the gravitational fit in Table 7, both nongravitational 
solutions based on the set {1,2,4, 5} improve the rnean 
residual only marginally (frorn ±1".79 to ±1".72/1".73), 
but reduce the offset from the reference apsidal line con- 
siderably (from 0°.078 to 0°.013/0°.006); (ii) both A\ and 
A 3 appear to be fairly well determined (with tlie relative 
errors of ±8-17%) and on the same order of magnitude, 
~10 -6 AU day -2 , as the parameters of tlie normal com- 
ponent of the nongravitational acceleration of the dwarf 
sungrazers C/2009 L5 and C/2006 J9 in Table 4 of Pa- 


a Dwarf Kreutz Sungrazer? 11 

per 1 ; and (iii) these parameters are one order of mag- 
nitude greater than the peak values of A\ for the cata- 
loged comets in nearly-parabolic orbits and at least two 
orders of magnitude greater than the typical values of 
A\ for such cornets (Marsden & Williams 2008). In sum- 
rnary, there is some — though, due to the limited data, 
not overwhelming — evidence that in terms of the non- 
gravitational effects in its orbital motion, C/1945 XI may 
sliare the properties of sorne dwarf Kreutz sungrazers. 

6.3. Modified Nongravitational Solutions 

The EXORB7 orbit-determination code employed by 
the second author allows one to vary arbitrarily all five 
parameters of the nongravitational law (9) — the expo- 
nents ?n, n, k, the scaling distance ro, and the normal- 
ization constant a. This option facilitates a more robust 
orbital experimentation, when the standard nongravita- 
tional model of Marsden et al. (1973) fails to provide sat- 
isfactory results. In Paper 1 we gained some experience 
with what we hereafter refer to as a modified nongravi- 
tational law g m od(u ?’o)j which is given by Equation (9), 
but while tlie three exponents remain unchanged from 
the standard law, the scaling distance varies broadly, al- 
ways entailing a change in a as well. 

The physics behind this sublimation-law experimenta- 
tion involves a parallelism between the scaling distance 
?’o and tlie snow line, a boundary of the zone in vacuum 
(or near vacuum), beyond which it is cold enough for a 
volatile substance to exist only in the solid phase. Cal- 
ibrated by water ice, for which in the isothermal model 
?’o ~ 2.8 AU, this distance depends in the first approxi- 
rnation on the effective latent heat of sublimation, L (in 


Table 9 

Orbital Solutions X/// {1, 2,4, 5} and {1, 2,4, 5} 
(Equinox J2000). 


Orbital element 


<{U,4,5} V S ^{1,2,4,5} 

Perihelion time t n (1945 TT) 

Dec. 27.97928 Dec. 27.97938 

Argument of perihelion uj 

67°.910 

67°. 883 

Longitude of ascending node Í2 

345°. 437 

345°.412 

Orbit inclination i 


141°.572 

141°.565 


/( AU) 

0.0071167 

0.0071250 


l (%) 

1.5287 

1.5307 

Orbit eccentricity 


0.99992638 

0.99992629 

Longitude of perihelion 

282°.82 

282°.83 

Latitude of perihelion B n 

+35°. 16 

+35°. 16 

Reference apsidal-line offset 

0°.0133 

0°.0060 

Parameter Ai (10 _5 AU day -2 ) 

+0.113 a 


Parameter A 3 (10 _5 AU day -2 ) 


+0.234 b 

Osculation epoch 1945 (TT) 

Dec. 11.0 

Dec. 11.0 

Mean residual (weighted) 

±1"73 

±1?72 

Distribution of Residuals O — 

C 

Time of 

W t ^{l,2,4,5} 

1-2, 4, 5} 

observation 







1945 (UT) 

R.A. 

Dec. 

R.A. Dec. 

Dec. 11.04687 

+4?9 

+2".0 

+479 +176 

12.04691 

-8.2 

-7.3 

1 

00 

o 

1 

to 

14.07000 

+0.3 

+ 1.0 

+0.2 +1.1 

15.06885 

-0.1 

-0.6 

0.0 -0.7 

a With a mean error of ±0.019 X 

10 -5 AU day -2 

(Table 8). 

D With a mean error of ±0.018 X 10 5 AUday 2 

(Table 8). 
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cal mol 1 ), of the sublimating species: 


const^ 2 

~J' 


( 10 ) 


where r 0 is expressed in AU and the constant is equal 
to 19100 AU"ï calmol -1 . Marsden et al. (1973) showed 
that in a plot of log (normalized sublimation rate) against 
log (heliocentric distance) the value of r 0 shifts the curve 
left or right along the axis of heliocentric distance r. 
Hence, by properly clioosing r 0 , the erosion-driven non- 
gravitational effects in the orbital rnotion of any cornet 
can approximately be expressed by a universal curve of 
log (normalized sublimation rate) against log(r/r 0 ). 

Highly refractory materials, such as metals or silicates, 
can sublimate only very close to the Sun, so that their 
snow line and scaling distance are rnuch smaller than 
1 AU. For example, for forsterite, the Mg end-member 
of the olivine solid solution system, tlie latent heat of 
sublimation is 130000 cal mol” 1 (Nagahara et al. 1994), 
so that its r 0 ~ 0.02 AU. 

In limiting cases the empirical equation (9) offers the 
following expressions for variations in the sublimation 
rates (equally applying to the modified and standard 
laws): 

lim g mod (r; r 0 ) ~ r -2 15 , 

x—ï0 

lim 5mod(r-;r 0 )~r -25 ' 65 . (11) 


where x — t/tq. The first lirnit approaches an extreme 
scenario in which all incident solar energy is spent on sub- 
limation, while the second lirnit crudely approximates the 
other extreme, when the energy is spent entirely on heat- 
ing the object (and increasing its thermal reradiation). 
The actual limiting expressions are r -2 for r/r 0 —> 0 and 


7mod(r-;r 0 ) ~ exp 



for r/r 0 —> oo, 


( 12 ) 


where the sublimation heat L is related to r 0 by Equation 

(10) , R is the gas constant, and T is the temperature at 
heliocentric distance r. The differences relative to (11) 
are due to the approximate nature of the law (9), which 
is responsible for a peculiar feature when applied to the 
orbital motion of C/1945 XI (see below). To illustrate 
it, we begin by introducing a local slope C of a modified 
nongravitational law between r and r + Ar; similarly to 

(11) , it is 


where 
C(r;r 0 ) = 


lim g mo d(r; r 0 ) ~ r 
Ar—ï0 


(13) 


dlnffmod(r;ro) nk 

dln(r/ro) m l + (r/r 0 ) -íl ’ 


The lirnits on </ are those given in (11). Since r 0 is nearly 
3 AU in the standard law, the contribution to the inte- 
grated nongravitational effect frorn r r 0 is always neg- 
ligible and the exact sliape of the standard nongravi- 
tational law at heliocentric distances r^>r 0 is unimpor- 
tant. However, when r 0 in a modified law is considerably 
smaller tlian tliat of the standard law, the power-law ap- 
proximation to the exponential law for the sublimation 



Figure 4. The nongravitational effects in the orbital motion of 
a comet derived from the modified law < 7 m od( r ; r 0 ) applicable to 
an arbitrary scaling distance ro. At heliocentric distances r > r* 
the law 7mod( r ! r o) matches the exponential variations in the sub- 
limation rate better than 5 , mo d( r ; r o)- The laws g mo d( r: > r o) an( i 
7mod( r ; r o) have equal logarithmic slopes at r^. The slope £ of 
<?mod( r ; r o) at r > r^ is practically constant. 


rate beyond ro could play a major role for comets whose 
perihelion distances q are much greater than ro. 

To find out at what heliocentric distance r^ does a 
modified law begin to deviate from the sublimation law 
7mod( r 5 r o) 5 we first define the local slope £ of 7 mo d( r 5 r o) 
similarly to (14): 


^ln7mod(r;r 0 ) _ Ly/r _ Ly/rp TF 
ain(r/ro) “ 2 RT 0 ~ 2 RT 0 \j r 0 ’ 


(15) 

where T(r) = T 0 /y/r at the heliocentric distances of ex- 
tremely low sublimation rates; for the isothermal nrodel, 
T 0 = 278.3 K. Since, from Equation (10), Ly/ró is a con- 
stant and L^/ró/RT 0 = 34.54, the point of deviation of 
the modified law ff m od(r;r 0 ) from the exponential sub- 
linration law, given by the condition of equal slopes, 
£(r*;r 0 ) =C( 7 +;r 0 ), cleternrines r # by comparing (14) 


with (15). After inserting the nunrerical values, 


?+ = 2.12 r 0 . 


(16) 
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Less than a millionth part of the Sun’s incident energy 
at the distance r is spent on the sublimation. 

Requiring that g mo d{r *; r 0 ) = Jm O d{r*-,r 0 ) and writ- 
ing the exponential sublimation law in the forrn 

7mod(r;r 0 ) = (3 exp 34.54 ^l - \jj-J > ( 17 ) 

we have 

/3 = 0.0267a. (18) 

The relationship between g m od(r;r 0 ) and 7 mod(r;r 0 ) is 
displayed in Figure 4. At present, the 7 mod(r;r 0 ) law is 
not incorporated into the orbit-determination code that 
we ernploy, as in rnost cases it would rnake hardly any 
difference numerically. However, the approximation of 
7mod(r; r 0 ) by g mo d(r; r 0 ) at r > r* (i) prevents one frorn 
testing nongravitational laws with a variable exponent / 
at these heliocentric distances and (ii) yields identical 
results from all solutions based on the scaling distances 
?’ 0 that are by rnore than 2.12 ?’ 0 smaller than the least 
heliocentric distance at which the cornet is observed, as 
seen frorn Figure 4. 

In Paper 1 we employed the modified-law paradigm to 
great advantage. We determined for all eight in-depth 
examined dwarf sungrazers that the astrometric posi- 
tions were fitted better (and the offsets from tlie ref- 
erence line of apsides came out to be smaller) when the 
nongravitational acceleration affecting their orbital rno- 
tions was assumed to vary with heliocentric distance r 
rnuch rnore steeply than prescribed by the standard law 
gstd{r)- The orbital motions of six out of the eight cases 
were fitted best with ?’ 0 <0.13 AU, which corresponds to 
L > 50 000 cal mol -1 , and for none of the eight was r 0 
higher than ~2.1 AU. The heliocentric distances of the 
dwarf Kreutz sungrazers in Paper 1 varied between 0.037 
ancl 0.068 AU (or between 8 and 15 Rq), so that these ob- 
jects were considerably closer to the Sun than C/1945 XI 
when under observation from Boyden. 

For C/1945 XI, the modified-law nongravitational so- 
lutions with ?’ 0 < 0.2 AU fit rnost satisfactorily the pre- 
ferred set of four observations, {3 0 } = {1,2,4, 5}. This 
rnatch is better than that from the standard-law solu- 
tions listed inTable 9. We refer to these modified-law so- 
lutions as A nlod [?’ 0 ]{3 ! 0 }, where X lras the same meaning 
as before and r 0 is in AU. Because the comet was not ob- 
served closer to the Sun than 0.6 AU, all tliese solutions 
are identical, as pointed out more generally under (ii) 
below Equation (18). We confirm this result by running 
solutions for r 0 equal to 0.01, 0.05, and 0.10 AU. We then 
select, as examples, the solutions AC od [0.05]{S 0 } and 
A/’ n ^ d [0-05]{S o } and present thern in Table 10; relative 
to the results based on the standard law, the offsets frorn 
the reference line of apsides now dropped, respectively, 
from 0°.013 to 0°.005 and from 0°.006 to 0°.003, whilc 
the weighted rnean residuals dropped, respectively, frorn 
±1"73 ancl ±1"72 down to ±1"66 in both R and N. 

The nongravitational parameters A\ and A 0 are now 
determined with relative errors of only ±5-7%. Their 
apparent discrepancy of nearly five orders of magnitude 
compared to those of the standard law is due entirely 
to the very different steepness of tlre two laws’ slopes. 
Since the comet was observed between 0.718 AU frorn 
the Sun (on December 11) and 0.598 AU frorn the Sun 


Table 10 


Orbital Solutions A//^ d [0.05]{ïío} an d A/’JJ od [0.05]{3 0 } 
(Equinox J2000). 

Orbital element A^ d [0.05]{3o } A^ d [0.05]{9o } 

Perihelion time t n (1945 TT) Dec. 27.97937 Dec. 27.97940 

Argument of perihelion u 67°.876 67°.869 

Longitude of ascending node Í2 345°.403 345°.397 

Orbit inclination i 141°,562 141°,561 

_ ., ^ f (AU) 0.0071256 0.0071279 

Perihelion distance q j (%) 1-5308 1. 5313 

Orbit eccentricity 0.99992629 0.99992626 

Longitude of perihelion L n 282°.83 282°.84 

Latitude of perihelion B n +35°. 16 +35°. 16 

Reference apsidal-line offset 0°.0055 0°.0034 

Parameter Ai (10 - 10 AUday -2 ) +0.307 a . 

Parameter A 3 (10 _ 10 AUday _2 ) . +0.877 b 

Osculation epoch 1945 (TT) Dec. 11.0 Dec. 11.0 

Mean residual (weighted) ±1"66 ±1?66 


Distribution of Residuals O 

-C 


Time of 

Af n ^ d [0-05]{9o} 

A^[0.05]{9b} 

1945 (UT) 

R.A. 

Dec. 

R.A. 

Dec. 

Dec. 11.04687 

+4?7 

+ 2?0 

+4'Í8 

+1"7 

12.04691 

- 8.0 

-6.9 

-7.9 

-6.9 

14.07000 

+0.3 

+0.9 

+0.3 

+ 1.0 

15.06885 

- 0.1 

-0.5 

0.0 

- 0.6 


a With a mean error of ±0.021 X 10 10 AU day 2 . 
b With a mean error of ±0.046 X 10 _10 AU day -2 . 


(on December 15), the magnitudes of the effective non- 
gravitational accelerations need to be compared in this 
range of heliocentric distances, rather than at 1 AU 
from the Sun, to test whether the results frorn tlre stan- 
dard law and the modified law are generally consistent. 
This test shows that the radial components from the 
two laws reach the same value, ±0.3 x 10 -5 AUday -2 , 
at r = 0.639 AU, while the norrnal components equally 
reach ±0.6 x 10 -5 AU day -2 at r = 0.648 AU, in either 
case well inside tlie observed range of heliocentric dis- 
tances. The comet’s favorable comparison with the dwarf 
sungrazers, mentioned in Section 6.2, is consistent with 
the results frorn the modified-law solutions: C/1945 XI 
appears to be subjected to a nongravitational accelera- 
tion tliat is near the lower end of the range of accelera- 
tions typical for the dwarf sungrazers. 

In summary, we obtain sorne evidence suggesting that 

(i) the nongravitational solutions fit the observations of 
C/1945 XI better than the gravitational solutions; and 

(ii) a law with a very steep slope (of the type ~r -25 ) may 
be superior to tlie standard law. Over all, however, the 
promising performance of the modified law has not over- 
turned our conclusion expressed at the end of Section 6.2 
that the limited data sample does not rnake a conclusion 
about the existence of large nongravitational effects in 
the rnotion of C/1945 XI sufficiently compelling. Every 
effort should be expended to improve the comet’s orbit 
by extending tlie observed arc with the help of additional 
images on Boyden patrol plates. 

6.4. Astrometric and Orbital Accuracy 

It may seem to be questionable to use the astrometric 
data of less than high quality for supporting the type of 
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Table 11 

Geocentric Ephemeris for Comet C/1945 XI from 1945 September 12 to 1946 January 30 (Equinox J2000). 


From solution X" I /)f d [0.05]{Sro} App. magnitude Differential ephemerides from solutions 

- Phase - - 


Date TT 


R.A. 

Dec. 

A 

r 

angle 

PA (RV) 

r- 4 

r -7 

^std 

{3o} 

S{3o} 

orbit D' 

1945 Sept. 

12 

7 h 39 m 0 

-io°oo' 

2.860 

2.466 

20° 1 

258°5 

15.7 

19.7 

s 

-0.31 

+ 2(5 

-2 S 8 

+117" 

s 

-0.25 

+10(3 


22 

7 50.1 

-11 48 

2.592 

2.310 

22.7 

263.5 

15.2 

19.0 

-0.31 

+2.7 

-2.8 

+116 

-0.26 

+10.3 

Oct. 

2 

8 01.6 

-14 03 

2.315 

2.148 

25.6 

268.2 

14.6 

18.2 

-0.32 

+2.8 

-2.9 

+115 

-0.27 

+10.4 


12 

8 13.7 

-16 54 

2.031 

1.980 

28.8 

272.8 

14.0 

17.3 

-0.33 

+3.0 

-3.0 

+113 

-0.29 

+10.4 


22 

8 27.0 

-20 36 

1.744 

1.804 

32.5 

277.2 

13.3 

16.3 

-0.34 

+3.2 

-3.1 

+110 

-0.31 

+10.4 

Nov. 

1 

8 42.3 

-25 33 

1.458 

1.619 

37.2 

281.5 

12.4 

15.1 

-0.35 

+3.4 

-3.2 

+107 

-0.35 

+10.3 


11 

9 02.4 

-32 30 

1.180 

1.423 

43.5 

285.4 

11.4 

13.6 

-0.38 

+3.5 

-3.3 

+94 

-0.40 

+10.0 


21 

9 35.1 

-42 48 

0.922 

1.213 

53.0 

286.9 

10.2 

11.9 

-0.42 

+3.2 

-3.5 

+71 

-0.53 

+8.9 

Dec. 

1 

10 54.6 

-57 49 

0.713 

0.982 

69.1 

277.1 

8.7 

9.7 

-0.46 

+1.5 

-3.3 

+30 

-0.90 

+5.6 


6 

12 34.0 

-65 04 

0.645 

0.855 

80.7 

256.5 

7.9 

8.5 

-0.36 

+0.1 

-2.1 

+8 

-1.22 

+1.7 


11 

15 11.8 

-65 27 

0.617 

0.719 

94.6 

219.8 

7.0 

7.0 

-0.03 

+0.1 

-0.1 

+1 

-0.97 

-3.4 


16 

17 12.4 

-56 11 

0.640 

0.569 

108.9 

189.8 

6.1 

5.3 

+0.04 

+1.7 

-0.3 

-4 

-0.39 

-5.9 

1946 Jan. 

5 

18 32.1 

-42 17 

0.662 

0.434 

126.4 

201.3 

5.0 

3.4 

+0.99 

-8.7 

+9.8 

+204 

+1.10 

+19.9 


10 

19 14.6 

-55 37 

0.566 

0.601 

114.8 

184.7 

6.1 

5.5 

+2.02 

-7.4 

+22.1 

+219 

+2.12 

+22.3 


20 

0 37.9 

-69 21 

0.525 

0.882 

84.7 

104.0 

7.6 

8.2 

+0.66 

+17.4 

+36.6 

+508 

+1.78 

+44.0 


30 

3 32.8 

-48 43 

0.664 

1.124 

60.6 

70.5 

9.1 

10.6 

-1.47 

+4.7 

-2.8 

+723 

-0.87 

+55.6 


in-depth orbital analysis of C/1945 XI tliat is described 
in the previous sections. With the observed arc of nrerely 
four days and the images taken with small cameras of 
poor spatial resolution and long exposures tracked at 
sidereal rates, the data may at first sight seem overex- 
ploited. However, the addressed objectives are not only 
to learn as much as possible about tlre orbital behavior 
of C/1945 XI in the context of other Kreutz sungrazers, 
but, as illustrated in Table 11, also to provide a basis for 
constraining an ephenreris for the comet’s nrotion both 
back and forward in tinre from tlre observed period, a ca- 
pability that is anrply put to use in our search for nrore 
images of the comet, as described in Section 7. 

7. SEARCH FOR MORE DATA ON COMET C/1945 XI 

The preceding sections demonstrate tliat, in spite of 
some evidence to the contrary, the current level of our 
knowledge of comet C/1945 XI is less than sufhcient for 
addressing the issue of whether or not tlris object was 
indeed a dwarf Kreutz sungrazer. It is essential that, 
where at all possible, existing infornration be rechecked, 
reanalyzed, and reinterpreted, and that a vigorous search 
for, and exanrination of, additional data be undertaken. 
All our efforts will be aided substantially by products of 
the Harvard College Observatory’s DASCH Project, cur- 
rently in progress, that we referred to in Section 3. As 
for the new data, our proposed work lras three objectives: 
(i) to identify, astronretrically measure, and reduce any 
further images of the conret in the collection of digitized 
Boyden patrol plates (once they beconre available) with 
the aim to refrne the orbit by extending the arc covered 
by observations as far back in tinre as possible; (ii) to 
constrain the conret’s physical behavior by deriving nrag- 
nitudes from all relevant digitized photograplrs and conr- 
pare the comet’s light curve with that of other sungraz- 
ers; and (iii) to search for and exanrine the conret’s post- 
perihelion inrages taken at Boyden for any traces of the 
object or its debris, including the presence, appearance, 
and orientation of a headless tail, diagnostic of the ob- 
ject’s fate near perihelion. 


7.1. Search for Additional Images 

Our extensive orbit exanrination in Section 6 allows 
us now to extrapolate, with a fair degree of confrdence, 
the comet’s motion way beyond tlre period of December 
11 15 and to realistically estimate uncertainties of an 
ephenreris as a function of tinre. The nominal geocentric 
ephenreris of tlre conret for tlre period of time fronr 1945 
Septenrber 12.0 TT to 1946 January 30.0 TT, presented 
in Table 11, is conrputed from the nongravitational solu- 
tion •A/’J^ 1 [0.05]{5o} (whose elenrents and residuals are in 
Table 10), but the differential ephenrerides, derived from 
the solution ■A/’jy^Qo} (elenrents and residuals in Ta- 
ble 9), fronr the gravitational solution {/{So} (Table 7), 
and from the orbit D' (see the text near Table 6), are 
also included for conrparison. The two nunrbers listed in 
Table 11 for each of these three additional solutions are 
the corrections in right ascension and declination, respec- 
tively, tlrat are to be added to the nonrinal coordinates 
in colunrns 2 and 3; note tlre different units used. 

The Hrst six colunrns of Table 11 are self-explanatory. 
The seventh column gives the position angle of the pro- 
longed radius vector (direction away fronr the Sun), and 
the eiglrth and ninth colunrns include a predicted ap- 
parent nragnitude on two assumptions for the brightness 
variations: as an inverse fourth and inverse seventh power 
of heliocentric distance, with no phase effect. Both ver- 
sions were normalized to the conret’s reported magnitude 
7 at discovery on Decenrber 11, which predicts a magni- 
tude of 9.5 at 1 AU from both the Sun and the Earth 
for the inverse fourth-power law and 10.6 for the inverse 
seventh-power law (see Section 7.2). 

We find the agreement among the ephemerides based 
on the orbit D' and the two nongravitational solutions 
cjuite remarkable. On the other hand, the ephenreris from 
the gravitational solution based on the sanre observations 
as the two nongravitational orbits is very disappointing, 
differing fronr the other tlrree by as much as 0°.2 over the 
period of 4.5 months. We note that the three consistent 
ephemerides conre fronr the orbital sets that imply that 
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Table 12 


Additional Boyden Plates with Possible Images of Comet C/1945 XI. 


Plate 

number a 

UT time at 
mid-exposure 

Plate center b 

Exp. 

time 

(min) 

Logbook 

reference 

Comet’s position (absolute and relative) b,c 

Observer d 

R.A. 

Dec. 

R.A. 

Dec. 

Plate 

Dist. 

P.A. 

Edge 

RB 14049 

1945 Sept. 19.11997 

8 h 04 m 6 

—15°17' 

60 

rblO.148 

7 h 46"b 

-11° 13' 

ON 

5891 

31286 

6885 

Bl(?) 

AM 25110 

Oct. 10.07683 

7 04.5 

-15 09 

85 

am45b_0144 

8 11.4 

-16 18 

ON 

16.11 

96.3 

1.18 

du Toit 

RB 14086 

17.07399 

8 04.0 

-30 17 

120 

rbl0_132 

8 20.3 

-18 39 

ON 

12.20 

18.5 

2.33 

j(?) 

AM 25132 

27.06607 

9 03.5 

-45 24 

60 

am45b_0148 

8 34.4 

-22 55 

OFF 

23.25 

342.8 

0.71 

Bl(?) 

AM 25145 

Nov. 5.05881 

9 04.2 

-30 24 

90 

am45b_0150 

8 49.7 

-28 04 

ON 

3.93 

305.4 

13.99 

Bester 

AM 25150 

12.06740 

10 04.0 

-45 29 

80 

am45b_0152 

9 05.1 

-33 25 

ON 

16.54 

311.8 

4.87 

j(?) 

RB 14121 

13.03559 

8 03.2 

-45 17 

120 

rbl0_138 

9 07.6 

-34 16 

OFF 

16.51 

53.8 

2.12 

Bl(?) 

AM 25169 

15.06371 

10 33.7 

-60 31 

85 

am45b_0154 

9 13.4 

-36 11 

OFF 

27.50 

323.1 

0.50 

J(?) 

AM 25187 

Dec. 7.04693 

11 04.6 

-45 32 

90 

am45b_0156 

13 04.0 

-65 58 

OFF 

25.94 

152.4 

1.49 

du Toit 

AM 25190 

8.01580 

10 33.7 

-60 31 

90 

am45b_0156 

13 34.2 

-66 28 

ON 

20.50 

126.1 

0.63 

j(?) 

AM 25208 

14.03336 

16 40.9 

-70 12 

40 

am45b_0158 

16 34.1 

-60 38 

ON 

9.60 

355.0 

11.94 

j(?) 

AM 25231® 

1946 Jan. 5.04600 

12 05.2 

-70 33 

90 

am45b_0162 

18 32.4 

- 42 24 

OFF 

52.65 

112.7 

31.37 

Bl(?) 

RB 14183® 

6.04985 

12 05.2 

-60 33 

115 

rblO.148 

18 38.1 

-44 57 

OFF 

55.57 

121.9 

35.99 

Bl(?) 

RB 14184 f 

8.07729 

17 17.5 

-47 19 

40 

rblO.148 

18 53.5 

-50 19 

OFF 

16.02 

109.7 

3.88 

Bester 

AM 25240 

21.79113 

0 05.1 

-44 27 

90 

am45b_0164 

1 36.3 

-66 22 

OFF 

25.13 

158.5 

1.89 

Britz 

AM 25241 

21.84446 

1 33.6 

-59 29 

60 

am45b_0164 

1 37.8 

-66 15 

ON 

6.78 

176.5 

14.73 

Britz 

AM 25244 

24.85635 

5 25.4 

-68 55 

90 

am45b_0164 

2 40.0 

-59 43 

OFF 

19.62 

277.3 

2.26 

Bl(?) 

AM 25260 

30.85589 

5 03.9 

-29 52 

90 

am45b_0166 

3 38.4 

-47 05 

ON 

23.86 

217.8 

2.56 

Bl(?) 


a AM plates taken with the Cooke lens (covering 34°48 in R.A. and 43°10 in Dec.); RB plates with the Ross-Fecker camera (22°32 by 27°90). 
b EquinoxJ2000. 

c Column Plate indicates whether or not the comet’s predicted position is within the plate limits; columns Dist. and P.A. show the predicted 
position relative to the plate center; and in column Edge is the distance of the comet’s image from the plate edge (either on or off the plate). 
b The identities of the observers whose abbreviations are B1 and J could not be determined, major efforts notwithstanding; see also Table 1. 
e This plate may show a segment of an early post-perihelion tail, if the comet still produced significant amounts of dust (see Figure 6 or 7). 

^ This is the only deliberate attempt at detecting C/1945 XI after perihelion; unfortunately, the searched position is incorrect (see Figure 5). 


the line of apsides lies within 0°.01 of the reference value, 
while the t/{3o} set leaves an offset of nearly 0°.08. 

To further appraise the magnitude of ephemeris uncer- 
tainties, we completed 200 Monte Carlo runs with pure 
gravitational solutions fitted to the {3o} set of observa- 
tions, assuming, conservatively, measuring errors of up 
to ±9" in the December 11-12 positions and up to ±2".5 
in the December 14-15 positions. The standard devia- 
tions from the 200 clones vary, in the 1945 September- 
November period, between 2'4 and 3' in right ascension 
and between 1'.4 and 5' in declination; in the 1946 Jan- 
uary period, the range is frorn 6' to nearly 20' in right 
ascension and from 8' to 31' in declination. On the av- 
erage the offsets of the t/{3o}-based ephemeris in Table 
11 amount to about 40% of these standard deviations. 
To sum it up, we believe that our nominal ephemeris is 
accurate, in a worst case scenario, to a few arcmin over 
the entire period of interest to us. 

The next step is a search for all patrol plates whose 
fields cover (or coulcl cover) the comet’s ephemeris posi- 
tion. Ideally, the pliysical size of each plate and its scale 
provide the angular distance frorn the center to an edge 
as a function of a position angle. It is in principle very 
easy to establish whether the comet’s image is or is not 
located within tlie plate’s limits. In reality, tlre problem 
is less straightforward because of a potentially significant 
error in tlre reported position of the plate center. Because 
of this uncertainty (much larger tlran in the eplremeris), 
the comet’s predicted position may not be within the lim- 
its of a plate when it should, or vice versa. The search 
is further complicated by the comet’s unknown bright- 
ness, the observing conditions (the sky transparency, the 
rnoon interference, the limiting magnitude, etc.), the ex- 


posure tirne, and the distribution of field stars (needed 
for the astrometry), all of which determines whether the 
comet’s irnage, even if present on the plate, can in fact be 
recognized, astrometrically measured and reduced, and 
photometrically evaluated. 

Table 12 presents information on a total of 18 plates 
that may show an irnage of C/1945 XI or its debris, 
including the unreduced plate AM 25208 frorn December 
14 (Table 1) but not tlre five from Table 4. Because of the 
position uncertainties, we include plates that nominally 
miss the comet’s position by up to 2°.5. 

These are all patrol plates, with very large fields of view 
(Section 3). The Ross-Fecker plates (RB) cover a rectan- 
gle of 22°.32 in right ascension and 27°.90 in declination 
(~600 square degrees), the Cooke plates (AM) 34°.48 by 
43°.10 (~1500 square degrees). The column of Table 12, 
titled Plate, shows whether the comet’s predicted posi- 
tion, given by the standard polar coordinates relative to 
the plate center, is (ON) or is not (OFF) within the limits 
of the plate. The colunm Edge indicates tlie distance of 
the comet’s predicted position from tlre plate edge. The 
smaller tlre number in an ON case, tlie more difhcult it 
will be to identify and astrometrically evaluate the im- 
age because of optical imperfections far frorn the optical 
axis of the instrument and a potential lack of appropriate 
field stars. The smaller tlie number in an OFF case, the 
greater is a cliance that the irnage could, after all, appear 
on the plate because of tlre positional uncertainties. 

We list a total of seven preperihelion plates on which 
the comet sliould show up if bright enough and four such 
plates with the comet’s predicted location barely rnissed; 
see Section 7.2 for tlre brightness constraints. After peri- 
helion tlrere are only two potentially useful plates by the 
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end of January to check whether any material part of the 
head survived; two additional plates just miss the comet’s 
expected position. On one of the plates taken on January 
21, we predict that tlie comet’s location is not far frorn 
the center. The comet’s head should also appear on five 
plates exposed during February (AM 25269 on the 7th, 
AM 25272 and AM 25273 on the 8th, RB 14224 on the 
19th, and RB 14246 on the 27th), but there is no point 
in examining these unless the comet shows up, ratlrer 
unexpectedly, on the January plates. 

Finally, there is a group of three plates (footnotes e and 
/ in Table 12), which rniss the predicted position of the 
comet’s head by a wide rnargin but may be relevant to a 
search for traces of the comet’s dust debris (Section 8.2). 
This group includes tlie January 8 plate (Table 1), taken 
witli the Ross-Fecker carnera specifically for the purpose 
to recover C/1945 XI (Section 3). Another plate, taken 
almost simultaneously with the Metcalf refractor (Ta- 
ble 1), is not listed, because it is deerned useless on ac- 
count of its small fielcl. 

7.2. Photometry and the Light Curve 

To predict the brightness of a cornet is always risky, as 
has amply been documented by historical examples. For 
the Kreutz sungrazers, the task is further complicated by 
apparent differences between tlre shapes of the preperihe- 
lion liglit curves of a prominent Kreutz cornet (observed 
as a bright object from the ground) and a dwarf sun- 
grazer (observed only frorn SOHO/STEREO ), as sug- 
gested by Kniglit et al. (2010). However, a preperihe- 
lion light curve is known for only one prominent sun- 
grazer, C/1965 S1 (Sekanina 2002). 8 It is not advisable 
to generalize this light curve as a standard attribute of 
all bright Kreutz sungrazers. This caveat is supported by 
the perceived differences arnong the post-perihelion light 
curves for five prominent Kreutz objects (C/1843 Dl, 
C/1882 Rl, C/1963 Rl, C/1965 Sl, and C/1970 Kl), 
which were fading with heliocentric distance at average 
rates from ?’ -3 - 3 to r -5,1 (Sekanina 2002). 

Knight et al. (2010) concluded that the rate of bright- 
ening of the dwarf Kreutz sungrazers observed witli the 
SOHO coronagraphs changes strikingly near 24 Rq (0.11 
AU) frorn the Sun frorn to ~ r -3 - 8 ; this slower 

rate of brightening then holds, according to them, down 
to 16 Rq (0.075 AU), where the light curve begins to 
flatten. They estimated tlrat tlre steep rate is unlikely to 
persist at heliocentric distances greater than 50 Rq (0.23 
AU) and that farther from tlie Sun the dwarf sungrazers 
brighten on their way to perihclion by again following 
an r -3 - 8 law. However, Ye et al. (2014) suggested that 
there is a significant diversity in the preperihelion light 
curves of the dwarf Kreutz sungrazers. They called at- 
tention to either an outburst (whose amplitude exceeded 
5 magnitudes) or an extremely steep briglrtening (r -11 
or steeper) of a dwarf sungrazer C/2012 E2 between 1.06 
AU and 0.52 AU frorn the Sun, while the object followed, 
on the average, an ~ r -4 (or flatter) law during its ap- 
proach to perihelion between at least 0.52 AU (rather 

8 Only a single, crude preperihelion magnitude estimate is avail- 
able for C/1882 Rl, the brightest Kreutz sungrazer over the past 
two centuries. This estimate is based on reports of its rivaling 
Venus in brightness at the time of discovery, 12 days before peri- 
helion (Gould 1883). 


than 0.23 AU) and 0.07 AU. Ye et al. also pointed out 
that another dwarf Kreutz sungrazer, C/2012 U3, was 
not detected at a heliocentric distance of 1.31 AU before 
perihelion and was then rnuch fainter than it should have 
been if it followed tlie Knight et al.’s (2010) prediction. 

Very illuminating is Ye et al.’s (2014) comparison of the 
light curves of C/2012 U3 and C/2011 W3 (Lovejoy). Al- 
though C/2011 W3 was not a dwarf sungrazer, it bright- 
ened according to an r~ 6 ’ 9 law (Sekanina & Chodas 2012) 
frorn 0.75 AU down to 0.34 AU, thus behaving very dif- 
ferently from the prominent sungrazer C/1965 S1 (r^ 41 
between 1.02 AU and 0.03 AU; Sekanina 2002). It ap- 
pears tlrat there is (i) a distinct possibility that, arnong 
the dwarf sungrazers, the steep, law applies to helio- 
centric distances of around 1 AU and possibly even fur- 
ther away frorn the Sun; and (ii) at least sorne of these 
objects appear to be subjected to outburst-like brighten- 
ing at moderate distances from the Sun on their way to 
perihelion. 

The potential implications for the brightness evolution 
of C/1945 XI are obvious. A rnajor obstacle to accept- 
ing tlris object as a dwarf sungrazer is its considerable 
brightness reported at the tirne of discovery. Compari- 
son with C/2011 W3 sliows that C/1945 XI was intrin- 
sically (i.e., after a correction lias been applied for the 
geocentric distance) fully 4.5 magnitudes (!) brighter at 
the sarne heliocentric distance. Assuming the reported 
estimate is correct, tlris problem can only be overcome 
if one accepts that C/1945 XI experienced a rnajor out- 
burst sorne tirne before the discovery, which was tliereby 
substantially facilitated. If the amplitude of the outburst 
was greater than ~4.5 magnitudes, and therefore cornpa- 
rable to tlrat of C/2012 E2, C/1945 XI would have been 
fainter than C/2011 W3 before the onset of the outburst. 
Frorn the standpoint of cliances that C/1945 XI could 
be detected on any of the pre-discovery plates in Table 
12, the pivotal parameter is the time when the putative 
event took place. If this general scenario is correct, the 
comet’s unusual brightness is indeed explained, but there 
is ratlier a slim cliance of detecting tlie comet on any of 
the pre-outburst plates. 

An alternative explanation of the reported magnitude 
of C/1945 XI, if it should be a dwarf sungrazer, is simply 
a gross overestimation of its briglrtness by the discoverer. 
In any case, tliese considerations show the critical need 
for a thorough photometric examination of the comet’s 
Boyden irnages, both those frorn the period of December 
11 15 and any pre-discovery ones. 

Keeping all our options open, we illustrate the dra- 
rnatic differences, over a wide range of heliocentric dis- 
tances, between the C/1945 XI brightness predictions 
based on the r -4 and r~‘ laws by listing the apparent 
magnitudes in colunms 8-9 of Table 11. One should, 
however, be able to discriminate between tlre two laws 
over a rnuch sliorter range of r. In fact, in the known 
Boyden irnages taken over the period of December 11- 
15, the comet should have brightened by 0.7 magnitude if 
the r~ 4 law applied, but by 1.3 magnitudes if tlre steeper 
law was in effect. The ephemeris suggests that C/1945 
XI miglrt have been as bright as magnitude 15 in tlie sec- 
oncl half of September, about 3 months before perihelion. 
In order that no irnage be rnissed, we began our search 
for plates that fit the predicted position of the cornet as 
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early as mid-September, since the Ross-Fecker carnera, 
the rnore powerful of the two patrol instruments, has 
nominally a limiting magnitude 15 as well (Section 3). 

8 . WHAT KIND OF A SUNGRAZER WAS C/1945 XI? 

Evidence of the fate of C/1945 XI after its passing 
through perihclion is vital for determining the place of 
this object in the hierarchy and evolution of the Kreutz 
system and for answering the question of whether it was 
a dwarf sungrazer. We propose to establish tliis frorn 
comparison with some better investigated sungrazers. 

8.1. C/1945 XI and the Transition Sungrazers 

Based on large numbers of observations, tlre members 
of the Kreutz system are usually divided into two cate- 
gories: prominent (bright and often quite spectacular as 
seen from the ground) and dwarf (defined in Section 1). 
This is a very simplistic classification, because it ignores 
intermediate objects, which occupy tlre transition be- 
tween both categories and, albeit scarce, are momentous 
for a better understanding of the disintegration process 
of the Kreutz sungrazers. 

In an effort to determine where C/1945 XI is likely to 
fit in, we selected three “standards” to define a scale for 
transition objects in the order of increasing similarity to 
the dwarf members. The standards were chosen to be 
represented by C/2011 W3, C/1887 Bl, and C/2007 L3. 
Next we describe their diagnostic properties, which will 
be searched for on the best suited plates of C/1945 XI 
in order to classify the cornet on tliis scale. 

C/2011 W3 was the most persistent and nearest the 
status of a prominent mernber among the three objects. 
With a perihelion distance of 1.19 Rq, tlris cornet sur- 
vived intact (without its nucleus disintegrating com- 
pletely) until 38 ± 5 hours after perihclion, at which tirne 
it suddenly fell apart in a terminal outburst, losing sud- 
denly its residual, ~10 12 g nucleus and head (Sekanina 
& Cliodas 2012). The tail, made up of microscopic dust 
particles ejected frorn the nucleus before perihelion, was 
seen to survive perihelion for at least 24 liours as a syn- 
dyname — a locus of dust released at different tirnes but 
subjected to the same radiation pressure acceleration /3 
of /3 = 0.6 the solar gravitational acceleration, typi- 
cal for dielectric submicron-sized particles. Only grains 
ejected less than 0.1 day before perihelion, which ap- 
proached tlre Sun within 1.8 Rq, sublimated away. After 
perihclion the cornet forrned a new, spectacular, head- 
less tail, rnost of which was a product of the event ~38 
hours after perihelion. This tail was a synchrone — a lo- 
cus of dust grains of different sizes (subjected to a broad 
range of radiation pressure accelerations), all released at 
the same tirne — that was under observation for up to 
90 days, until mid-March 2012. Its representative length 
during rnost reported sightings again suggested a maxi- 
rnurn radiation pressure acceleration of /3 ~ 0.6 the solar 
gravitational acceleration. Only during tlie late obser- 
vations (60-90 days after perihelion) did the visible tail 
consist of merely micron-sized and larger grains, whose 
/3 <C 0.6 the solar gravitational acceleration. 

C/1887 B1 is the runner-up to C/2011 W3. Avail- 
able information is ratlrer limited, because tliis object 
was discovered only after perihelion as a headless, nar- 
row, ribbon-like tail, again a synchrone. Nonetlieless, the 
systematic variations in this tail’s orientation with tirne 
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suggest that the nucleus survived intact until 5.8 ± 0.8 
hours after perihelion (Sekanina 1984), at whicli point it 
rnust have suddenly disintegrated in a terminal outburst 
similar to tliat of C/2011 W3. At perihelion the cornet 
was merely 1.04 Rq from the Sun (Sekanina 1978; Seka- 
nina & Chodas 2004). The tail, truly spectacular during 
the early sightings from 8.5 days after perihclion on, was 
under observation until 18.5 days after perihelion and its 
length implied a maximum radiation pressure accelera- 
tion of dust particles that dropped frorn ~0.4 the solar 
gravitational acceleration at discovery to ~0.05 during 
the last three reported observations. The minimum par- 
ticle size in the visible part of the tail was thus increasing 
with tirne frorn less than 1 rnicron to several microns. 

C/2007 L3 is, of tlie three standards, the rnost re- 
mote from the prominent Kreutz-system members, be- 
ing in fact a representative of a relative srnall group of 
bright dwarf Kreutz sungrazers, which develop long, ex- 
tremely narrow tails upon their approacli to perihclion 
and which we refer to hereafter as the superdwarfs. Like 
ordinary dwarf sungrazers, their nuclei fail to survive per- 
ihelion, so that they forrn no post-perihclion tails. We 
chose C/2007 L3 as a representative primarily because it 
was imaged by both SOHO and STEREO, unlike most 
of its peers. The comet had a perihelion distance of 
1.530 Rq (Marsden 2008), very close to our best esti- 
mate for C/1945 XI, its brightness peaked at magnitude 
~3 and its head disappeared about 80 minutes before 
perihelion, at a heliocentric distance of 2.75 Rq (Green 
2007). The SOHO and STEREO images were studied 
in considerable detail by Thompson (2009), who closely 
conhrmed the earlier conclusions by Sekanina (2000), 
based on several objects, that preperihelion images of 
these tails show that they are synchrones, referring to 
emission-shutoff times of 26 ± 7 hours before perihelion 
and to heliocentric distances at shutoff of 23.5 ± 4.5 Rq. 9 
This average agrees remarkably well with Knight et al.’s 
(2010) distance at which the slope of the light curves sud- 
denly changes frorn —7.3 to —3.8 (Section 7.2). 10 The 
tail of C/2007 L3 survived the comet’s perihelion, as did 
the preperihelion tail of C/2011 W3. Thompson exam- 
ined its remnant and found, unlike Sekanina & Chodas 
(2012) for C/2011 W3, that the feature remained syn- 
chronic to its last examined image, nearly 17 hours af- 
ter perihelion, but that its position was slightly mod- 
ihed by a loss of angular momentum, possibly due to 
atmospheric drag from the solar corona. The appar- 
ent discrepancy is explained by the fact that the dust 
emission of C/2011 W3 did not shut off about 1 clay 
before perihelion, but continued. Both Sekanina (2000) 
and Thompson (2009) detected the maximum radiation 
pressure accelerations of /3 ~ 0.6 the solar gravitational 
attraction. Comet C/2007 L3 was not the only Kreutz 
superdwarf whose dust tail was observed to survive per- 
ihelion. Other examples are C/1979 Ql, the hrst Sol- 
wind comet (Michels et al. 1982), C/1998 K10, one of 
the dozen comets examined by Sekanina (2000), etc. It is 

9 For the best two imaged comets (C/1996 Y1 and C/1998 K10) 
in his set, Sekanina (2000) noticed that the tail was actually bent 
and, unlike its bright synchronic portion, the faint, outer part was 
a syndyname with a range of the ejection times. 

10 This evidence may be interpreted to indicate that dust emis- 
sion essentially ceases at 24 Rq from the Sun and the r~ 3-8 law of 
brightening measures an effect of the object’s fragmentation. 
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likely that most, if not all, superdwarfs would show their 
tails surviving perihelion, if their irnages are closely ana- 
lyzed. About 10 such cornets from 1995-2005 in Knight 
et al.’s (2010) Table 2 and all 19 sungrazers from 2006- 
2013 in Sekanina & Kracht’s (2013) Table 1 belong to 
this subcategory, whiclr is at present estimated to include 
about tlrree dozen objects. 

So where does C/1945 XI fit in? Employing the con- 
straint that no tail was ever reported with the unaided 
eye after perihelion, it is certain that this comet was too 
faint to be placed between C/2011 W3 and C/1887 B1 
or even near C/1887 Bl. Since a preperihelion tail is 
always found to survive for only 1 day or so after peri- 
helion, there is no cliance of detecting it on any of the 
plates listed in Table 12. Sliould we find traces of the 
comet’s synchronic tail, a product of a post-perihelion 
outburst, on any of these plates, then C/1945 XI is to 
be classified on this three-point scale between C/1887 B1 
and C/2007 L3, in which case it was a transition object 
more massive than a superdwarf. If not, it has to be 
concluded that C/1945 XI cannot be positioned higher 
than C/2007 L3, so tliat, at best, it was a superdwarf. 

8.2. Search for Traces of C/191^5 XI After Perihelion 

From these considerations it follows that while a search 
for, and examination of, the comet’s preperihelion im- 
ages on the plates listed in Table 12 should serve to im- 
prove our understanding of the comet’s orbital rnotion 
and its light curve, a careful search for traces of the cornet 
on, and examination of, the plates taken after perihelion 
should provide information on the comet’s fate and help 
answer the question of whether C/1945 XI was a transi- 
tion object or a dwarf sungrazer; because of tlre lack of 
information on tlie comet’s appearance in close proximity 
of perihelion (Table 12), we cannot distinguish between 
its classification as a dwarf or a superdwarf. 

Before we discuss specific searcli opportunities, we note 
that the dust released frorn tliis comet before perihclion 
was affected by sublimation less severely than the dust 
frorn C/2011 W3. The orbital elements for C/1945 Xl’s 
preperihelion ejecta subjected to radiation pressure ac- 
celerations /3 = 0.6 the solar gravitational acceleration 
are summarized in Table 13. It follows that if this 
comet’s dust had sublimation properties similar to those 
of the dust from C/2011 W3, only particles ejected within 
~1 hour of perihclion should have sublimated completely, 
because of a greater perihclion distance of C/1945 XI. 
A major delay of several days in passing through perihe- 
lion is noted as a yet another effect: for example, parti- 
cles ejected 50 days before the comet’s perihelion did not 
reach tlieir own perihelion until late on January 5. 

As to the sources of information on comet C/1945 XI 
after perihelion, a rernote opportunity occurred on 1946 
January 3.51 UT, the tirne of a solar eclipse. 11 Unfor- 
tunately, it was merely partial, of magnitude 0.553, and 
visible only frorn Antarctica and surrounding oceanic re- 
gions. Even less of the Sun was occulted at the locations 
of several permanent British scientific stations in Antarc- 
tica near the southernmost tip of Soutli America (Oper- 
ation Tabarin; e.g., Haddelsey 2014). At J.Shanklin’s 
suggestion we checked the UK Met Office Archives in 

11 See the website http://eclipse. gsfc.nasa.gov/SEplotl901/ 
SE1946Jan03P.GIF. 


Table 13 

Orbital Elements of Dust Particles Subjected to /3 = 0.6 
and Ejected from Comet 1945 XI up to Perihelion. 


Time of 
ejection 
(days) a 

Time of 
perihelion 
(days) a 

Perihelion 

distance 

(%) 

Orbital 

eccentricity 

-50 

+8.868 

3.757 

1.035 

-40 

+7.142 

3.749 

1.041 

-30 

+5.410 

3.732 

1.049 

-25 

+4.540 

3.721 

1.055 

-20 

+3.667 

3.706 

1.064 

-15 

+2.788 

3.682 

1.077 

-10 

+1.899 

3.641 

1.101 

-7 

+1.357 

3.596 

1.127 

-5 

+0.989 

3.545 

1.157 

-3 

+0.610 

3.450 

1.217 

-2 

+0.412 

3.356 

1.279 

-1 

+0.203 

3.134 

1.426 

-0.7 

+0.136 

3.027 

1.527 

-0.5 

+0.089 

2.896 

1.642 

-0.3 

+0.041 

2.673 

1.861 

-0.2 

+0.017 

2.481 

2.082 

-0.1 

-0.004 

2.144 

2.567 

-0.05 

-0.010 

1.843 

3.150 

-l h 

-0.009 

1.779 

3.299 

0.0 

0.000 

1.530 

4.000 


a Reckoned from the time of the comet’s perihelion passage. 


Exeter and, not surprisingly, found no report on any rel- 
evant observations at the time. 

Back at Boyden, an attempt was made by M. J. Bester 
on January 8 to recover the comet; the search, as pointed 
out (Sections 3 and 7.1), was conducted with the Metcalf 
refractor and the Ross-Fecker piggyback camera. Both 
exposures pointed in the same direction because they 
overlapped each other; the Metcalf 45 minute exposure 
began 20 minutes earlier and was longer by 5 minutes. 
The pointing of the telescope is, however, a mystery: as 
performed, it could not accommodate the comet’s pre- 
dicted position within the limits of either plate. 

January 8 was the date of publication of Cunningham’s 
(1946b) ephemeris, but early in the morning Bester could 
not liave been in the possession of the official announce- 
ment from Copenhagen. Yet, the telescope was set to 
point precisely in accord with the epliemeris in decli- 
nation and was — most peculiarly — off by exactly 
2 hours to the west of the ephemeris place in right as- 
cension. Sucli a coincidence surely cannot be accidental. 
We presume that Bester was notified directly by Cun- 
ningham aliead of the ephemeris’ publication, perhaps 
by cable, but there was a misconnnunication, an error 
in one coordinate. In any event, the time of Bester’s 
observation was fundamentally incorrect because during 
the exposures tlie comet, less than 30° from the Sun, 
was below and, later, very close to the horizon. Since 
the comet was almost exactly south of the Sun, it could 
have been observed either shortly after sunset or, some- 
what preferably, just before sunrise. The need to take 
a fairly long exposure did of course constrain the choice 
of tlie observation time. If observed at the correct time 
witli the telescope set to point in accord with Cunning- 
ham’s ephemeris in both coordinates, the comet’s po- 
sition would have been exposed on both plates, even 
though only ~0°.5 from the edge of tlie Metcalf plate. 
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ANGULAR DISTANCE FROM COMET IN RIGHT ASCENSION 
Figure 5. The predicted projected position of C/1945 XI relative to the sky area exposed on the plate RB 14814, taken by M. J. Bester 
with the Ross-Fecker camera on 1946 January 8.07729 UT. While the comet’s head was more than 4° off the eastern edge of the plate, 
much of the dust tail made up of particles released near and shortly before perihelion occupy the lower, southern part of the plate’s field. 
Dust ejected as far back as ~40 days before perihelion occupies the rest of the plate. The curves to the south of the plate center are 
the synchrones, the loci of the dust that was ejected from the comet at particular times between 1.2 days before perihelion and 1.6 days 
after perihelion; the curve is the perihelion synchrone. The numbers along the synchrones are the values /3, the ratio between the solar 
radiation pressure acceleration and the solar gravitational acceleration; /3 and the ejection time govern the motions of the dust ejecta. 
Values near /3 ~ 0.1 are typical for micron-sized grains, those near or in excess of 0.5 refer to a variety of submicron-sized grains. The thick 
curve at the top, labeled /3 = 0.6 is a syndyname defined by this /3 ratio and typical for dielectric submicron-sized particles. Shown along 
this curve are the times of ejection in days before perihelion; see text for the significance of this syndyname. 


Figure 5 displays the projected positions of the comet 
and the Sun relative to the pointing direction of the plate 
RB 14184, taken with the Ross-Fecker camera on Jan- 
uary 8, 11 days after perihelion. Although dust ejecta are 
predicted to be all over the plate, their detection, contin- 
gent on the time-dependent emission rate, is problematic 
on account of their generally low surface brightness, ex- 
cept when released during a sharply peaked event, such 
as a terminal outburst. The curves emanating in the 
figure from the comet’s head and distributed mostly in 
the lower part of the plate are the near-perihelion syn- 
chrones. Their selected range is from 1.2 days before per- 
ihelion to 1.6 days after perihelion. Few post-perihelion 
ejecta project onto the plate, barely touching its south- 
eastern corner, and only a limited fraction of the perihe- 
lion ejecta crosses the plate. On the other hand, rnuch 
of the dust released as far back as ~40 days before per- 
ihelion is projected onto the plate. We emphasized in 
Section 8.1 that the preperihelion tails of the transition 
object C/2011 W3 and of C/2007 L3 and many other 
(if not all) superdwarfs survived perihelion for not much 
longer than ~1 day at the most. Such features are too 


close to the Sun to observe from the grouncl. It is only for 
the sake of interest that Figure 5 shows the syndyname 
/3 = 0.6, rnade up of dielectric submicron-sized particles 
ejected from C/1945 XI before perihelion, to miss the 
field of the sky covered by the plate RB 14184, nearly 
grazing it along the entire northern edge. 

The plate RB 14184 may prove useful. If C/1945 XI 
experienced a terminal outburst and disintegrated within 
an hour or two after perihelion and if the ejecta contained 
a large abundance of grains several microns in diameter, 
there is a possibility that this debris could show up as 
a band of light across the lower left corner of the plate. 
The phase angles favor a modest forward-scattering ef- 
fect, being mostly in a range of 100°-110°. Submicron- 
sized ejecta released from the comet at the same times 
should miss the plate. Thus, a detection or nondetection 
of a post-perihelion synchronic tail could provide rnean- 
ingful constraints on the existence of a terminal outburst 
and on the tirne of nucleus disintegration. 

We found two rnore plates that might provide addi- 
tional limited information on a potential post-perihelion 
terminal outburst: AM 25231, taken with the Cooke lense 
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BOYDEN EXPOSURE AM 25231 
TAKEN ON 1946 JANUARY 5 
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Figure 6. Predicted positions, on a Cooke plate AM 25231 taken 
on January 5.04600 UT, of synchrones that depict C/1945 Xl’s 
dust ejecta from hypothetical terminal outbursts 2, 3.5, and 6 hours 
after perihelion. The smaller-font digits refer to particles subjected 
to a variety of radiation pressure accelerations /3 between 0.34 and 
1.05 the solar gravitational acceleration. Since the Kreutz comets 
appear to have /3 < 0.6, only outbursts that occurred between ~2 
and ~4 hours after perihelion might be detected on this plate. 

in the morning of January 5, ancl RB 14183, taken with 
the Ross-Fecker camera 24 hours later. They both are 
listed in Table 12. The synclirones for the January 5 
plate are depicted in Figure 6, tliose for the January 6 
plate in Figure 7. The positions of dust particles rela- 
tive to the comet’s nucleus, determined by a computer 
code in terms of their polar coordinates, were converted 
to the polar coordinates of the offsets frorn the plate cen- 
ter, using a technique described in some detail in tlie Ap- 
pendix. An advantage of these plates relative to that on 
RB 14184 is their exposure times, 115 and 90 minutes, re- 
spectively, as opposed to only 40 minutes (Table 12). The 
January 5 plate might provide useful information only if 
the terminal outburst took place between, at most, ~2 
and ~4 hours after perihelion; a synclrronic tail contain- 
ing earlier ejecta would miss this plate, while segments of 
later synchrones limited to /3 > 0.6 the solar gravitational 
acceleration are unlikely to refer to any real tails (Sec- 
tion 8.1). Similarly, the January 6 plate rnight provide 
information only on an outburst that occurred between 
~1 and ~5 hours after perihelion. No forward-scattering 
effect can be expected in either case, as the phase angles 
range between about 60° and 80°. 

In the highly unlikely case tlrat a major fragment of 
the original nucleus survived for weeks after perihelion, 
the best search opportunity is offered by a Cooke plate 
AM 25241, taken on January 21, on which the comet is 
predicted to be less than 7° from the center (Table 12). 
Less favorable circumstances should accompany a search 
on a plate AM 25260, taken 9 days later, on whiclr the 
comet’s predicted position is only 2°.5 frorn the edge. Ad- 
ditional search opportunities are rnore discouraging still. 


BOYDEN EXPOSURE RB 14183 
TAKEN ON 1946 JANUARY 6 
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Figure 7. Predicted positions, on a Ross-Fecker plate RB 14183 
taken on January 6.04985 UT, of synchrones depicting the comet’s 
dust ejecta from hypothetical terminal outbursts 2, 3, and 5 hours 
after perihelion. The smaller-font digits refer to particles subjected 
to a variety of radiation pressure accelerations (3 between 0.3 and 
0.8 the solar gravitational acceleration. Since the Kreutz comets 
appear to have /3 < 0.6, only outbursts that occurred between ~1 
and ~5 hours after perihelion might be detected on this plate. 

9. CONCLUSIONS 

Comet C/1945 XI has been rather an oddball among 
the Kreutz sungrazers observed from tlre ground ever 
since its discovery 70 years ago. Although the dilatory 
handling of the Boyden plates contributed to the snail 
pace in getting out the facts about the object, it was 
nonetheless its apparent lack of luster near and after per- 
ihclion that is primarily responsible for our ignorance of 
its properties, fate, and place in the Kreutz system. 

The comet’s disappointing post-perihclion performance 

- especially the absence of a prominent dust tail pro- 
vides a strong argument against its being on a par with 
the headless sungrazers C/2011 W3 and C/1887 B1 that 
survived perihelion witlr a fairly massive nucleus intact 
but disintegrated shortly afterwards. The plausible sce- 
narios for C/1945 XI are thus limited to two: either the 
comet’s surviving rnass was sufficient to generate a mod- 
est, but not spectacular, post-perihelion tail, or the cornet 
completely disintegrated still before perihelion. While 
there is not much of a difference between the two scenar- 
ios, we use the second one to define a dwarf Kreutz sun- 
grazer. The existence of a terminal outburst is therefore 
critical to distinguish between tlre two options. 

Taking account of tlie indirect planetary perturbations 
on the premise that a precursor of C/1945 XI separated 
frorn its common parent with C/1882 R1 (Sekanina & 
Chodas 2004), we determine that the line of apsides of 
C/1945 XI should point toward an ecliptical longitude 
of L n = 282°.44 and latitude of B n = +35°. 16 (Equinox 
J2000) and that the comet’s osculating semimajor axis a 
should be 96.70 AU at an epocli 1945 December 28.0 TT 
or 96.67 AU at an epocli December 11.0 TT. 
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Employing this value of a as an orbital constraint, 
our Best Fit purely gravitational solution based on all 
five observations left an unacceptably large offset of 
0°.57 from the nominal apsidal-line direction, apparently 
because the third observation was inferior. This was 
confirmed by four-observation gravitational solutions, of 
which the one based on the positions from December 11, 
12, 14, and 15 was by far the rnost promising, leaving 
an offset from the nominal apsidal-line direction of 0°.08, 
much better but still not entirely satisfactory. 

Applying a technique introduced recently by us in Pa- 
per 1, we show that the incorporation of a nongravita- 
tional acceleration of the orbital motion improves the fit 
in terms of both the apsidal-line offset and rnean resid- 
ual. Use of Marsden et al.’s (1973) standard formalism 
suggests an acceleration on the order of 10 -6 AU day -2 
(with a relative error of ±8-17%), near the lower end 
of a range typical for the dwarf sungrazers. A prefer- 
able nongravitational model, based on a rnuch steeper 
acceleration law (apparently comrnon arnong the dwarf 
sungrazers; see Paper 1) than is the standard law, leads 
to the same magnitude of the integrated effect but fits 
the observations better, with the relative errors of the 
nongravitational parameters reduced to ±5-7%. 

In order to augment the current, very limited data set, 
an extensive search for furtlrer possible images of the 
comet should be undertaken. We hope to be able to 
conduct such a search, based on our orbital computa- 
tions suggesting that at least 7 and perhaps as rnany 
as 11 preperihelion Boyden patrol plates taken between 
late September and mid-December may contain suclr im- 
ages. The plates will be inspected in the near future, once 
their digitized copies are made available by the DASCH 
Project. At present we do not anticipate to expand our 
search to archives of wide-field plates from other obser- 
vatories (e.g., Tsvetkov & Tsvetkova 2012). 

In addition to positional data of the cornet, its plio- 
tometry will also have to be performed to reexamine the 
reported brightness at discovery and to settle the issues 
of the rate of brightness variations witli heliocentric dis- 
tance and tlre cliance of an early preperihelion outburst. 
Comparison with C/2011 W3 and C/2012 E2, arnong 
others, may hclp settle some of the issues of preperihe- 
lion activity of the Kreutz sungrazers. 

Finally, to pursue a solution to the complex problem 
of C/1945 XI’s place in the hierarchy of the Kreutz 
system and the comet’s fate, we also will examine some 
post-perihelion Boyden patrol plates to search for both 
potential relics of the cornet itself and/or possible traces 
of its dust tail as a product of a modest terminal out- 
burst that may have occurred just hours past perihelion. 
Wliereas it appears that, masswise, C/1945 XI could not 
compete with either C/2011 W3 or C/1887 Bl, we will 
try to disentangle a mystery of whether it still coulcl be 
considered a transition object distinctly rnore massive 
than a superdwarf, such as C/2007 L3, C/1998 K10, 
or C/1979 Ql, or it does not differ materially from the 
Kreutz dwarfs or superdwarfs — and tlius answer the 
cardinal question posed in this paper’s title. 
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APPENDIX 

TRANSFORMATION OF DUST PARTICLE POSITION FROM 
COMETOCENTRIC TO PLATE CENTER COORDINATES 

The code computing the position of a dust particle, which 
is a function of the ejection time, í e ject, and a ratio of 
the radiation pressure acceleration to the solar gravita- 
tional acceleration, /3, provides the polar coordinates in 
reference to the comet’s nucleus in projection onto the 
plane of the sky — the angular distance from the nu- 
cleus, D* (í e je C t, /3), and the position angle, H^ (í ejec t, /?), 
reckoned from the north tlirough the east. The equatorial 
coordinates of the nucleus, right ascension a 0 and decli- 
nation do, and the coordinates of a particle, a* (í e ject, /3) 
and 5*(t e j ec t,0), are related to and n* by the well- 
known expressions, 

cos D* = sin Sq sin ± cos Jo cos 6 % cos Aa t , 

cot n* = cos Jo tan 5* csc Aa^. — sin Jo cot Aa*, (19) 

where Aa* = a* — ao- To convert the particle’s position 
from the cometocentric coordinate system to the coordi- 
nate system centered on the center of a plate, we first 
compute right ascension a* and declination ó* of the 
particle from the formulas 

sin 2 D + tan n* (sin ± cotD* cos Jo sec n*) 

* cos 2 D* (1± cos 2 J 0 tan 2 n^.)— sin 2 <3 0 

( 20 ) 

and 

tan6 * = cosAa^.tam5o(l±cotn^tanAa^. cscJ 0 ). (21) 

Since the sign of cos Aa* lias a direct effect on the sign 
of tanJ^. and therefore also on the signs of sinJ^. and 
cos D *, the quadrant of Aa* in Equation (20) needs to 
be chosen such that after inserting the values of Aa* and 
6 * from (21) into the hrst equation of (19) one gets the 
correct value of D* and not its supplement. 

Next, identifying a 0 and J 0 with the coordinates of 
the plate center, rather than the comet’s nucleus, and 
inserting them together with the particle’s coordinates 
a* and ó* into Equations (19), one obtains an angular 
distance D* and a position angle n* that determine the 
particle’s position relative to the plate center. 
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